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ABSTRACT 

We present analysis of the near- infrared spectra of 114 rest-frame UV-selected star-forming galaxies 
at z '^ 2. By combining the Ha spectra with photometric measurements from observed 0.3-8 ^m, 
we assess the relationships between kinematics, dynamical masses, inferred gas fractions, and stellar 
masses and ages. The Ha line widths give a mean velocity dispersion a — 101 ± 5 km s~^ and a 
mean dynamical mass Mdyn ~ 6.9 ± 0.6 x 10^" M0 within a typical radius of ^ 6 kpc, after excluding 
AGN. The average dynamical mass is ~ 2 times larger than the average stellar mass, and the two are 
correlated at the 2.5a level and agree to within a factor of several for most objects, consistent with 
observational and systematic uncertainties. However, ~ 15% of the sample has Mdyn ^ M^. These 
objects are best fit by young stellar populations and tend to have high Ha equivalent widths Whq ^ 
200 A, suggesting that they are young starbursts with large gas masses. Rest-frame optical luminosity 
and velocity dispersion are correlated with 4ct significance; the correlation and its accompanying large 
scatter suggest that the processes which produce the correlation between luminosity and velocity 
dispersion in local galaxies are underway at 2; ~ 2. Fourteen of the 114 galaxies in the sample 
have spatially resolved and tilted Ha emission lines indicative of velocity shear. It is not yet clear 
whether the shear indicates merging or rotation, but if the galaxies are rotating disks and follow 
relations between velocity dispersion and circular velocity similar to those seen in local galaxies, our 
observations underestimate the circular velocities by an average factor of ~ 2 and the sample has 
(Vc) ~ 190 km s~^ . Using the local empirical correlation between star formation rate per unit area 
and gas surface density, we estimate the mass of the gas associated with star formation, and find 
a mean gas fraction of ~ 50% and a strong decrease in gas fraction with increasing stellar mass. 
The masses of gas and stars combined are considerably better correlated with the dynamical masses 
than are the stellar masses alone, and agree to within a factor of three for 85% of the sample. The 
combination of kinematic measurements, estimates of gas masses, and stellar population properties 
suggest that the factor of ~ 500 range in stellar mass across the sample cannot be fully explained by 
intrinsic differences in the total masses of the galaxies, which vary by a factor of ~ 40; the remaining 
variation is due to the evolution of the stellar population and the conversion of gas into stars. 
Subject headings: galaxies: evolution — galaxies: high redshift ~ galaxies: kinematics and dynamics 
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1. INTRODUCTION 

A galaxy's mass is one of its most fundamental prop- 
erties. Both popular models of galaxy formation and 
recent observations suggest that mass is an important 
factor in determining when and how quickly galaxies 
form their stars; star formation in more massive galaxies 
is b oth expected and observed to start at earlier times 



Dave et al.lf2005HHeavens et al.ll200l iJuneau etlll 



I2005|r The redshift range 1.5 < ^ < 2.5 has been 
shown to be particularly important both for the buildup 
of stellar mass and for accretion onto massive black 
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holes. A large fraction of the stellar mass in the uni- 
verse today likely formed at z > 1 (jDickinson et aQ 
(20031 iRudnicket al]l200"3(l . and this re dshift range also 
sees the peak of brigh t QSO activity l|Fan et alJl20flll: 
iDi Matteo et al.ll2003D . Effective techniques now exist 
for the selection of galaxies at z ^ 2; these use the galax- 
ies' observed optical ( Stcidcl ct al. 2004) or near-infrared 



(Daddi et al. 



LI 2003) colors, or a combination of the two 



2004!|1 , and can be used to select both star- 



forming and passively evolving galaxies. 

Even with large galaxy samples mass is difhcult to mea- 
sure, especially at high redshift. While galaxy formation 
models most naturally parameterize galaxies as a func- 
tion of mass, observers can only determine mass through 
the measurement of luminosity or (preferably but with 
more difficulty) kinematic properties. Measurements of 
galaxy mass at high (and low) redshift take a variety 
of approaches. The most direct methods use kinematic 
tracers to measure the depth of a galaxy's potential well. 
At relatively low redshifts, when rotation curves can be 
traced to large radii and galaxies' mass distributions 
modeled, this technique yields relatively robust measure- 
ments, as shown by the good agreement between re- 
sults derived fr om optical emission lines and 21 cm map- 
ning of H T re.g.lCourtea,un997HKobulnickv fc Gebh^TdH 
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I2OOOI) . At high redshift the situation is more comph- 
cated. While spatially resolved "rotation curves" can 
sometimes b e obtained usi ng nebular emission lines in 
the near-IR (iPettini et all [2 001: Erb et al. 2003, 2004| 
iForster Schreiber et al.l 120061) . the limitations imposed 
by the seeing mean that they do not necessarily pro- 
vide a unique kinematic model. A more robust measure- 
ment is the velocity dispersion, which can be measured 
for most high redshift galaxies and is less affected by 
the seeing. Velocity dispersions have been used to esti- 
mate the dynamical masses of relatively small samples of 
galaxies at 2 ~ 2 and z ~ 3; results typically range from 
~ iq i " to ~ 10" M^ llPettini et a,1 .■2001^ ' Erb et a,1pM 
12004 Ivan Dokkum et alJl2004 lSwinbank et al. 2003). A 
model for the mass distribution is still required to turn 
a measurement of the velocity dispersion into a mass es- 
timate, and this is a significant source of uncertainty in 
such dynamical masses. 

Dynamical masses measured from velocity dispersions 
do not, of course, trace the full halo masses of high red- 
shift galaxies; detectable nebular line emission is typi- 
cally confined to sites of active star formation within a 
central region of a few kpc. Halo masses can be esti- 
mated less directly, however, through the galaxies' clus- 
tering properties and the predicted clustering properties 
of halos as a function of mass. Adelberger et al. (2005b) 
find that the correlation lengths of the 2; ~ 2 rest-frame 
UV-selected galaxies on which this paper focuses corre- 
spond to halos with a typical mass of ^ 10^^ Mq. 

Estimations of stellar mass at high redshift are in- 
creasingly common, driven by wide-field IR detectors 
which allow large samples of rest-frame optical photome- 
try, and by the Spitzer Space Telescope which enables 
detection of rest-frame near-IR light. Stellar masses 
of z ~ 2-3 galaxies arc determined through popula- 
tion synthesis models applied to such broadband pho- 
tometry (e.g. Qp apovich et al. 2001; Shaolev ct al. i2001t 
IForster Schreiber et al.-2004:.Shaplcv et al.-2005. among 
many others); the masses obtained range from ~ 10^ to 
> 10"'^^ M0. Such masses usually represent lower limits 
to the true stellar mass, however, since it is possible for 
the light from older populations of stars to be obscured 
by current star formation. 

The most difficult component of a galaxy's baryonic 
mass to measure at high redshift is its gas. For the most 
luminous, gas-rich, or gravitationally magnified exam- 
ples, this can be estimated with millimete r observations 
(e.g. IRaker enni2004t iGreve et a,1.ll2fiol . but for the 
typical 2 ~ 2 galaxy such observations are still out of 
reach. As the gas probably accounts for a significant 
fraction of the baryonic mass in young galaxies at high 
redshift, we have made use of the empirical correlation 
between star formation rate density and gas surface den- 
sity to estimate the gas masses of the galaxies considered 
here. 

This paper is one of several making use of a large sam- 
ple of Ha spectra of 114 2 ~ 2 galaxies. In this work 
we discuss the galaxies' kinematic properties^ and th eir 
stellar and dynamical masses, while lErb et alJ l)2006bf) is 
devoted to star formation. This paper is organized as 
follows. We describe the selection of our sample, the ob- 
servations, and our data reduction procedures in |21 In 
ij^lwe outline the modeling procedure by which we deter- 
mine stellar masses and other stellar population param- 



eters. 21 is devoted to the galaxies' dynamical masses 
as derived from the Ha line widths; we compare them 
with stellar masses in M.l\ and assess the relationship 
between velocity dispersion and rest-frame optical lumi- 
nosity in t l4.2l We compare the distributions of velocity 
dispersions of galaxies at 2 ~ 2 and 2 ~ 3 in t l4.3| In 
i|Slwe examine galaxies with spatially resolved and tilted 
emission lines. [El describes our estimates of gas masses 
from the local correlation between star formation rate 
and gas density per unit area, and comparisons of stellar, 
gas and dynamical masses. We summarize our conclu- 
sions and discuss our results in fJJl Separately, we use 
the same sample of Ha spectra to construct composite 
spectra according to stellar mass to show that there is a 
strong correlation between increasing oxygen abundance 
as me asured by the [N II] /Ha ratio and increasing stellar 
mass IjErb e t al. 20063) • Galactic outflows in this sample 
are discussed by Steidel et al (2006, in preparation). 

A cosmology with Hq = 70 km s^^ Mpc~ , Qm — 0.3, 
and r^A = 0.7 is assumed throughout. In such a cosmol- 
ogy, 1" at 2 = 2.24 (the mean redshift of the current 
sample) corresponds to 8.2 kpc, and at this redshift the 
universe is 2.9 Gyr old, or 21% of its present age. For 
calcu lations of stellar masses and star formation rates we 
use a lChabrierl (2003) initial mass function (IMF), which 
results in stellar masses a nd SFRs 1.8 tim es smaller than 
would be obtained with a lSaloeted l)1955|) IMF. 

2. SAMPLE SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

The galaxies discussed herein are drawn from the rest- 
fra me UV-selected 2^ 2 spectroscopic sample described 
by iSteidel et al.l 1)2004(1 . The candidate galaxies were 
selected by their UnGTZ colors (from deep optical im- 
ages discussed bv ISteidel et al.ll2004|) . with redshifts for 
most of the objects in the current sample then con- 
firmed in the rest-frame UV using the LRIS-B spectro- 
graph on the Keck I telescope. The total sample of 114 
Ha spectra consists of 75 new ob servations and 39 that 
have been published previously l|Erb et alJ 12003", "200^ 
Shaolev et al. 2004). In total we attempted 132 galax- 
ies with previously known redshifts, but 23 were not de- 
tected. We also attempted to observe 28 galaxies with- 
out previously known redshifts but which met the pho- 
tometric selection criteria; only 5 of these yielded se- 
cure Ha redshifts. Most of the undetected galaxies have 
K > 21.5, while most of those we detect have K < 21.5; 
we are less likely to detect Ha emission for objects that 
are faint in K. This is not surprising, given the corre- 
lation between K magnitude and star formation rate re- 
cently found for the 2 '^ 2 sample by Rcddv et al. (2005iL 
This co rrelation is confirmed and discussed bv iErb et alJ 
i2006bH 

The set of galaxies presented here is not necessarily 
representative of the UV-selected sample as a whole, be- 
cause objects were chosen for near-IR spectroscopy for a 
wide variety of reasons. Criteria for selection included: 

1) galaxies near the hue of sight to a QSO, for studies of 
correlations between galaxies a nd metal systems seen i n 
absorption in the QSO spectra l|Adelberger et al.l2005a|) : 

2) morphologies — elongated in most cases, with a few 
more compact objects for comparison; most of these have 
been previously discussed by Er b et al.i.2004: 3) galaxies 
with red or bright near-IR colors or magnitudes, or oc- 
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casionally those whose photometry suggested an unusual 
spectral energy distribution (SED; 13 obje cts selected be- 
cause they have K < 20 are discussed bv lShaplev et al.l 
|2004); 4) galaxies that have excellent deep rest-frame 
UV spectra, or are bright in the rest-frame UV; 5) in 
order to confirm classification as AGN; or 6) members of 
close pairs with redshifts known to be favorable relative 
to night sky lines. The remaining objects are galaxies 
within ~ 25" of our primary targets, which were placed 
along the slit and observed simultaneously with the pri- 
mary targets. A significant number of galaxies meet more 
than one of these criteria. The sightline and secondary 
pair objects are therefore generally representative of the 
UV-selected sample as a whole, while the addition of 
the IR red or bright objects means that the total Ha 
sample is somewhat biased toward more massive galax- 
ies (see ^. We identify AGN as galaxies that show 
broad and/or high ionization emission lines in their rest- 
frame UV spectra, or as objects with broad Ha lines or 
very high [N II]/Ha ratios. The AGN fraction (5/114) 
is similar to the fraction found for the full spe ctroscopic 
sample in the GOODS-N field bv iReddv et' al. (2005), 
using direct detections in the 2-Ms Chandra Deep Field 
North images. 

In order to assess the representativeness of the NIR- 
SPEC sample, in Figure H we compare the colors and 
magnitudes of the galaxies in the NIRSPEG sample with 
that of all spectroscopically confirmed UV-selected galax- 
ies in the same fields and redshift range. The top panels 
show 7^ vs. G — TZ (upper left) and G — TZvs. Un — G (up- 
per right). The shaded regions in the upper right panel 
show the selection criteria for, from top to bottom , the 
z ~ 3 C/D and MD galaxies (iSteidel et alJl2003 D and 
the z ~ 2 BX and BM galaxies l|Steidel et alJl2004D . Ob- 
jects that fall outside the selection windows were either 
selected using earlier photometry or criteria extending 
slightly blueward in Un — G. The lower panels show Ks 
vs. TZ — Ks at lower left, and Kg vs. J — Kg at lower 
right. These comparisons show that the NIRSPEG sam- 
ple spans nearly the full range in TZ and Kg magnitudes 
and in TZ — Kg and J — Kg colors. We have observed 
~ 30-50% of the galaxies with NIRSPEG at the bright 
and red ends of the distributions, and ^ 10% at the faint 
and blue ends. Most objects selected by the BX/BM cri- 
teria have G — 7<^ < 0.5, as do all but three galaxies in 
the NIRSPEG sample; the reddest and bluest galaxies in 
G — TZ are somewhat underrepresented. 

For the purposes of comparisons with other surveys, 
we note that 10 of the 87 galaxies for which we have Ha 
spectra and J Kg photometry have J — Kg > 2.3 (the se - 
lection criterion for the FIRES survev. lFranx et all2003|) : 
this is similar to the ^12% of UV-selected galaxies which 
meet this criterion (jReddv et al. ,2005i) . 18 of the 93 
galaxies for which we have K magnitudes have Kg < 20, 
the selection criterion for the K20 survey IjGimatti et al.l 
12002') : this is a higher fraction than is found in the full 
UV-selected sample (~ 10%), bec ause we int entionally 
targeted many ii'-bright galaxies IjShaolev et al. 20041 . 
Four of the 10 galaxies with J — Kg > 2.3 also have 
Kg < 20. 

2.1. Near- 1 R Spectra 

All of the Ha spectra (with the exceptions of 
GDFb-BN88, Q0201-B13, and SSA22a-MD41, which 



were observed with the ISAAG spectrogra ph on the 
VLT and previously discussed bv lErb et al.l 120 03) were 
obtained using the near-IR spectrograph NIRSPEG 
(jMcLean et alJll998(l on the Keck II telescope. Observ- 
ing runs were in May 2002, May, July and September 
2003, and June and September 2004. The conditions 
were generally photometric with good (~ 0.5-0.6") see- 
ing, with the exceptions of the May 2003 and June 2004 
runs which suffered from occasional clouds and seeing up 
to ~ 1". Ha falls in the iiT-band for the redshift range 
2-0 S ^ S 2.6, which includes the vast majority of the 
current sample; 6 of the 114 objects have 1.4 ^ z < 1.8, 
and were observed in the if-band. We use the 0.76" x 42" 
slit, which, in NIRSPEG's low resolution mode, provides 
a resolution of - 15 A (- 200 km s"^ i? ~ 1400) in the 
iiT-band. 

Our observing procedure is described by lErb et all 
()2003(1 . For a typical object, we use four 15 minute inte- 
grations, although the number of exposure varied from 2 
to 6 for total integrations of 0.5 to 1.5 hours. We perform 
blind offsets from a nearby bright star, returning to the 
offset star between each integration on the science target 
to recenter and dither along the slit. In most cases we 
attempt to observe two galaxies with separation < 25" 
simultaneously by placing them both on the slit; thus the 
position angle is set by the positions of the two galaxies. 

The spectra wer e reduced using the standard proce- 
dures described bv lErb et alJ l)2003j) . Uncertainties were 
accounted for via the creation of a two-dimensional frame 
of the statistical la error appropriate to each pixel. The 
last step in the reduction was to extract one-dimensional 
spectra of each galaxy; this was done by summing the 
pixels containing a signal along the slit. The same aper- 
ture was then used to extract a variance spectrum from 
the square of the error image described above; the square 
root of this is a Icr error spectrum, which was used to de- 
termine the uncertainties in the line fluxes and widths. 
Line centroids, fluxes and widths were measured by fit- 
ting a Gaussian profile to the line using the IRAF task 
splot, with the base fixed to the average value of the 
continuum (we do not detect significant continuum flux 
for most of the objects in the sample). The splot task 
provides la statistical uncertainties for each of the fitted 
quantities, which are determined by a series of Monte 
Garlo simulations which perturb the spectrum according 
to the la uncertainties determined from the error spec- 
trum described above. The objects observed are listed in 
Tabled and some representative examples of t he spectra 
are sh own i n Figured Others ar e presented bv lErb et alJ 
ll200l and lShaplevetaLl ll200l . 

2.2. Near-IR Imaging 

J-band and Kg-hand images were obtained with the 
Wide-field IR Gamera fWIRG: Iwilson et al.ll2?X)l on 
the 5-m Palomar Hale telescope, which uses a Rockwell 
HgGdTe Hawaii-2 2kx2k array, with a field of view of 
8.5' X 8.5' and spatial sampling of 0.249 arcseconds per 
pixel. Observations were conducted in June and Octo- 
ber 2003 and April, May and August 2004. Some images 
were also taken in June 2004, courtesy of A. Blain and J. 
Bird. We used 120-second integrations (four 30-second 
coadds in K, one 120-second coadd in J), typically in 
a randomized ~ 8" dither pattern of 27 exposures. To- 
tal exposure times and 3a image depths in each field are 
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Fig. 1. — Photometric comparisons of the NIRSPEC sample with all spectroscopically confirmed C/nG7?.-selected galaxies in the same 
fields and redshift range. The top panels compare the rest-frame UV photometry: TZ vs. G — 7?. at upper left, a nd G — TZ vs. Un — G at 
upper right. From top to bottom, the shaded regions show selection criteria for the z ~ 3 LBGs (C/D and MD. ISteidel et ani2003f) and 
for BX and BM galaxies iSteidel et al. 2004). Objects that fall outside the selection windows were either selected using earlier photometry 
or criteria extending slightly blueward in Un — G. In the bottom panels we compare the near-IR photometry: Ks vs. 1Z — Kg at lower left, 
and Ks vs. J — Ks at lower right. The dashed lines show Ks = 20 and J — Ks = 2.3 in the left and right panels respectively. We include 
objects with 1.4 < z < 1.9; 6 of the 114 objects in the NIRSPEC sample are in this redshift range. 



given in Table [21 

Only those images with seeing as good as or better 
than our optical images (typically ~ 1.2"; 0.85" in the 
Q1700 field) were incorporated into the final mosaics. 
The images from each dither sequence were reduced, 
registered and stacked using IDL scripts customized for 
WIRC data by K. Bundy (private communication). The 
resulting images were then combined and registered to 
our optical images using IRAF tasks. Flux calibration 
was performed with reference to the 15-30 2MASS stars 
in each field. We have also used a 45 min image of 
Q2346-BX404 and Q2346-BX405 from the near-IR cam- 



era NIRC (iMatthews &: SoifeillT99^ on the Keck I tele- 
scope taken on 1 July 2004 (UT), courtesy of D. Kaplan 
and S. Kulkarni; this image was reduced similarly, and 
flux-calibrated by matching the magnitudes of objects 
in the field to their magnitudes in the calibrated Q2346 
WIRC image. 

Photometry was performed as described by 
iShaplev et al.l 1)2005(1 . Briefly, galaxies were de- 
tected (i.e., found to consist of connecting pixels with 
flux at least three times the sky a and area greater 
than the seeing disk) in the 7^-ba i id ca talogs as de- 
scribed in detail by ISteidel et alJ l(2003j) . with total 







rHE STELLAR, GAS AND DYNAMICAL MASSES OF ■. 


; - 2 GALAXIES 




40 




:HDF-BX1439; 


60 

40 


1 '; 


:Q2346-BX482. 


80 
60 




:Q1623-BX502: 


100 
80 
60 


r ': 


: 1 1 1 1 1 1 _ 
:Q1700-BX490: 


20 



1 ! 1 


k^ 


m 


20 



M 

I ,: 


fJH^ 


i^ 


40 

20 




y 


ftvrViV^ 


40 

20 




l(d 


UM 







2.08 


2.12 






30 


1 ': 


r 1 1 1 1 . , 

: Q2343-MD59- 


20 


- ': 






- 


10 


- : 


jl 


1 1 




- 





rt 


nJi 


*i 


i# 


i A 


"|l| \:\ '■ 1 


TT;,- 



2,24 2.28 1.96 2 

Wavelength (/zm) 



2.04 




2.24 



Fig. 2. — Representative examples of the Ha spectra. Each row shows objects drawn from each quartile of Ha flux: the top row contains 
galaxies drawn from the highest quartile in flux, the second row from the second highest quartile, the third row from the third highest 
quartile, and the bottom row from the lowest quartile. The dotted lines in each panel show, from left to right, the locations of [N IIJA6548, 
Ha, [N II]A6583, [S II]A6716, and [S II]A6732. 



TZ magnitudes determined by increasing the isophotal 
detection aperture by a factor of two in area. TZ — K and 
TZ — J colors were determined by applying the 7?.-band 
detection isophotes to the J and K images; the J and 
K magnitudes were then determined from the total TZ 
magnitudes and the TZ — K and TZ — J colors. We also 
applied the detection algorithm to the /-iT-band images, 
and then applied the resulting isophotes to the J-band 
images. This procedure resulted in two sets of J and 
K magnitudes which typically differed by ~ 0.3 mag 
or less, comparable to the photometric uncertainties 
derived from the simulations described below; these 
differences arise because the 7^-band total aperture may 
not include all of the K-hand light, or vice versa. For 
the final J and K photometry, we adopted the more 
significant of the two magnitudes for each galaxy. The 
additional detecti on in K represent a c hange in method 
from that used bv lShaolev et alJ l)2005j) : this change, as 
well as a slight change in the photometric zeropoints 
from new flux calibrations, means that some magni- 



tudes given here may be slightly d i fferen t than those 
previously published (|ShaDlev et alJ 120041 HiM) . We 
have also trimmed noisy, less well-exposed edge regions 
of the i^-band images somewhat more than previous 
versions for cleaner detections in K^ with the result 
that some objects with previously published magnitudes 
are no longer in the K-hand sample. Photometric 
uncertainties were determined by adding a large number 
of fake galaxies of known magnitudes to the images, and 
detecting them in both TZ and K bands to mimic our 
photometry of actual galax i es. Th is process is described 
in detail bv IShaolev et alJ l|2005|) . As noted in g] the 
photometric apertures used here are approximately the 
same size as those used to determine dynamical masses 
from the Ha spectroscopy. 

2.3. Mid-IR Imaging 

Two of our fields have been imaged by the In- 
frared Array Camera (IRAC) on the Spitzer Space Tele- 
scope. Images of the Q1700 field were obtained in Oc- 



6 



ERB ET AL. 



tober 2003 during the "In-Orbit Checkout", and have 
been previously disc ussed in detail (|Barmbv et alJl2()0l: 
iShaplev et al.ll2005|) . We also make use of the fully re- 
duced IRAC mosaics of the GOODS-N field, which were 
made public in the first data release of the GOODS 
Legacy project (M. Dickinso n, PI). Thes e images are 
described in more detail by iReddv et al.1 (j^05). For 
both fields we us e the photometric procedure described 
bv Shaplc v"erall (|2005). 

3. MODEL SEDS AND STELLAR MASSES 
3.1. Modeling procedure 

We determine best-fit model SEDs and stellar popu- 
lation parameters for the 93 galaxies for which we have 
iiT-band magnitudes. Most of these (87) also have J-band 
magnitudes, and 35 (in the GOODS-N and Q1700 fields) 
have been observed with the IRAC camera on the Spitzer 
Space Telescope. For modeling purposes we correct the 
K magnitudes for Ha emission; the typical correction 
is 0.1 mag, but for 4/93 objects (Q2343-BX418, Q2343- 
BM133, Q1623-BX455, Q1623-BX502), it is > 0.4 mag. 
We use a modeling procedure id entical to that described 
in detail bv lShaplev et alJ l)2005|) . and review the method 
briefly here. 

Photome try is fit w ith the solar metallicity 
iBruzual & Chariot! l)2003|) models; as shown by 
lErb et al.. (2006a), solar metallicity is a reasonable ap- 
proximation for the massive galaxies in the sample, and 
more typical galaxies have metallicities only slightly less 
t han solar . Although we use the models which employ 
a ISalpeteil lfl95 5) IMF, it is well known that the steep 
faint-end slope of this IMF below 1 Mq overpredicts the 
mass-to- light ratio (M/L) and stellar mass by a factor 
of ~ 2 (e.g. iBell et al.ll200l lRenzinill2005l) . For more 
accurate comparisons of stellar and dynamical masses 
we have therefore con verted all stellar masses and star 
formation rates to the iChabrien l)2003() IMF by dividing 
by 1.8 (a change in the faint end of the IMF affects only 
the inferred total stellar mass and star formation rate, 
because low mass stars do not contribute significantly to 
the flux in the observed passbands). Note that previous 
model SEDs of some of the same objects discusse d here 
used t he Salpeter IM F llShaplev et al.ll200ll2005|) . 

The iCalzetti et al.l l)2000j) starburst attenuation law 
is used to account for dust extinction. We employ a 
variety of simple star formation histories of the form 
SFR oc e(-*=f/^), with r = 10, 20, 50, 100, 200, 500, 
1000, 2000 and 5000 Myr, as well as r = oo (i.e. con- 
stant star formation, CSF). We also briefly consider more 
complex two component models, as described below; but 
as it is generally difficult to constrain the star formation 
histories even with simple models, we use only the single 
component models for the main analysis. For each galaxy 
we consider a grid of models with ages ranging from 1 
Myr to the age of the Universe at the redshift of the 
galaxy, with extinctions ranging from E{B — V) — 0.0 
to E{B — V) = 0.7, and with each of the values of r 
listed above. We compare the UnGTZJK+IRAC magni- 
tudes of each model (shifted to the redshift of the galaxy) 
with the observed photometry, and determine the nor- 
malization which minimizes the value of x^ with respect 
to the observed photometry. Thus the reddening and 
age are determined by the model, and the star forma- 



tion rate and total stellar mass by the normalization. In 
most cases all or most values of r give acceptable val- 
ues of x^i while the constant star formation models give 
the best agree ment with star for mation rates from other 
indicators fsee lErb et al.""2006b"l. For these reasons we 
adopt the CSF model unless it is a significantly poorer 
fit than the r models. CSF models are used for 74 of 
the 93 objects modeled, and r models for the remain- 
ing 19; as discussed by Shaplev et al. (2005), the most 
massive galaxies, with ^I^, > 10^^ M0, are usually sig- 
nifica ntly better fit by d eclining models with r = 1 or 2 
Gyr. lErb et all l|2006b() discuss this issue further. The 
adopted parameters and values of r may be found in Ta- 
ble El 

Parameters derived from this type of modeling are 
subject to substantial degeneracies and systematic un- 
certainties, primarily because of the difficulty in con- 
straining the star formation histo ry. These issues have 
been discussed in detail els ewhere l|Papovich et alJl2001l: 
IShaplev et "all 120011 120051). Here we simply note that 
age, reddening, and the value of r suffer most from these 
degeneracies, and that the stellar mass is less sensitive 
to the assumed star formation history and thus more 
tightly constrained. We examine the effects of the as- 
sumed star formation history, and of photometric errors, 
through a series of Monte Carlo simulations which we use 
to determine uncertainties on the fitted pa rameters. The 
simula tions are conducted as described bv lShaplev et alJ 
IJ2005D : the observed colors of the galaxy are perturbed 
by an amount randomly drawn from a Gaussian distri- 
bution with width determined by the photometric errors, 
and the best-fit SED is determined as usual for the per- 
turbed colors, including the star formation history r as 
a free parameter. We conduct 10,000 trials per object to 
estimate the uncertainties in each fitted parameter, and 
iteratively determine the mean and standard deviation 
of each parameter for each object, using 3(T rejection to 
suppress the effects of outliers. The resulting mean frac- 
tional uncertainties are {ax/ {x)) = 0.7, 0.5, 0.6 and 0.4 
in E{B — V), age, SFR and stellar mass respectively. The 
distributions of ax/{x) are shown in Figure |21 The addi- 
tion of the mid-IR IRAC data significantly improves our 
determination of stellar masses; for those objects with 
IRAC data (all those in the Q1700 and HDF fields) we 
find (Tm^Z-M^, ^ 0.2, while those without the IRAC data 
have aMt/Mi, ~ 0.5. The IRAC data also reduce the 
uncertainties in t he other param eters, though not by as 
large an amount. IShaplev et al.l (^2005-') show that while 
the IRAC data significantly reduce the uncertainties in 
stellar mass, the lack of these data does not substantially 
change or bias the inferred mass itself, especially when 
the K-h&nd magnitude can be corrected for Ha emission. 

Confidence interval plots for M* vs. E{B — V) and age 
vs. E{B — V) for selected objects are shown in Figure^ 
where the dark contours represent 68% confidence levels 
and the light 95%. The red x marks our adopted best 
fit. The objects shown are representative of the range of 
properties spanned by the NIRSPEC sample. The first 
two objects, Q1623-BX502 and Q2343-BX493, have low 
stellar masses, young ages and relatively high values of 
E{B - V). The third galaxy, Q1700-BX536, has prop- 
erties close to the average of the sample, and also shows 
how much smaller the confidence intervals can become 
with the addition of the IRAC data. The last object, 
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Q2343-BX610, is a massive, old galaxy which can only 
be fit with Ah > 10" M© and age tsi > 1 Gyr. 
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Fig. 3. — The distribution of uncertainties of each parameter 
fitted by our SED modeUng code, from the Monte Carlo simula- 
tions described in the text. We show histograms of the ratio of the 
standard deviation to the mean of each parameter for all the ob- 
jects modeled, determined iteratively with 3cr rejection. From left 
to right and top to bottom, we show E{B — V), age, star forma- 
tion rate, and stellar mass. Stellar masses are the best-constrained 
parameter, followed by the age of the star formation episode. 



3.2. Model SED results 

The best-fit parameters from the SED modeling are 
given in Table |21 and their distributions shown in the 
histograms in Figure |S| The mean stellar mass is 3.6 ± 
0.4 X 10^° Mq, and the median is 1.9 ± 0.4 x lO^" Mq. 
The mean age is 1046 ± 103 Myr, and the median age 
is 570 ± 137 Myr. The sample has a mean E{B - V) 
of 0.16 ± 0.01 and a median of 0.15 ± 0.01. The mean 
SFR is 52 ± 10 Mq yr-\ while the median is 23 ± 3 
Mq yr~^; the difference between the two reflects the fact 
that a few objects are best fit with high SFRs (> 300 
Mq yr~^). We discuss the res ults of the models further 
in the following sections and in lErb et alJ lj2006bjl . where 
we compare them with properties determined from the 
Ha spectra. 

The 93 galaxies whose Ha spectra and SEDs are dis- 
cussed here are a subset of a larger sample of 461 objects 
with model SEDs. This larger sample will be discussed 
in full elsewhere (Erb et al. 2006, in preparation); for the 
moment we compare the distribution of stellar masses 
of the NIRSPEC galaxies with that of the full sample, 
which is representative of the UV-selected sample except 
that it excludes the ^ 20% of objects which are not de- 
tected to K ^ 22.5 (as discussed by Erb et al. (2006b), 
these are likely to be objects with low stellar masses and 
relatively low star formation rates). For the purposes of 
comparison with the full sample, we use models for the 
NIRSPEC galaxies in which we have not corrected the 
K magnitude for Ha emission (since we do not have Ha 
fluxes for the full sample); everywhere else in this paper, 
we use corrected K magnitudes for the modeling as de- 
scribed above. We also use CSF models for all galaxies 



for this comparison. The distributions of stellar mass for 
the two samples are shown in Figure The mean stel- 
lar mass of the subsample for which we have Ha spectra 
is 3.7 ± 0.5 X 10^° Mq; this is shghtly higher than that 
of the full UV-selected sample which has a mean stellar 
mass of 3.0 ± 0.5 x 10^° Mq, though the two are consis- 
tent within the uncertainties. A slightly higher average 
mass is expected for the NIRSPEC sample, since some 
of the galaxies were selected because of their bright K 
magnitudes or red TZ — K colors. 

3.3. Two component models and maximum stellar 
masses 

One limitation of this type of modeling is that it 
is sensitive to only the current episode of star forma- 
tion; it is possible for the light from an older, un- 
derlying burst to be completely obs c ured by c urrent 
star formation fPapovi ch et al.ri200H lYan et all l2004t 
IShaplev et al. 2005). Although we can only weakly con- 
strain the star formation histories of individual galax- 
ies even with simple declining models, more complex 
two-component models are a useful tool to assess how 
much mass could plausibly be missing. Two types of 
two-component models were used to examine this ques- 
tion. The first is a maximal mass model, in which we 
first fit the X-band (and IRAC, when it exists) data 
with a nearly instantaneous (t = 10 Myr), maximally 
old (zform — 1000) burst, subtract this model from the 
observed data, and fit a young model to the (primarily 
UV) residuals. Such models produce stellar masses ^ 3 
times higher than the single component models on aver- 
age, but for galaxies with (single component) M+ < 10^° 
Mq, the two component total masses are 10-30 times 
larger. Such models are poorer fits to the data, how- 
ever, and are not plausible on average because the young 
components require extreme star formation rates; the av- 
erage SFR for the young component of these models is 
~ 900 Mq yr"^, > 30 times higher than the average star 
for mation rate predi cted by the Ha and UV emission 
fsee lErb et al.l l2006b') . Nevertheless, we cannot rule out 
the possibility that such models may be correct for some 
small fraction of the sample. 

We also use a more general two-component model, in 
which we fit a maximally old burst and a current star 
formation episode (of any age and with any of our val- 
ues of r) simultaneously, allowing the total mass to come 
from any linear combination of the two models. As ex- 
pected given the number of free parameters, such models 
generally fit the data as well as or better than the single 
component models. They do not, however, usually signif- 
icantly increase the total mass; the total two component 
mass is more than three times larger than the single com- 
ponent mass in only six cases. In general, then, the data 
do not favor models with large amounts of mass hidden 
in old bursts. 

4. VELOCITY DISPERSIONS AND DYNAMICAL MASSES 

The Ha velocity dispersion a reflects the dynamics of 
the gas in the galaxies' potential wells. Because it re- 
quires only a measurement of the line width, it can be 
determined for most of the objects in the sample; it is 
therefore our most useful kinematic quantity. Out of 
the 114 galaxies with Ha detections, we measure a ve- 
locity dispersion from the Ha line width for 85 objects. 
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Fig. 4. — Sample confidence intervals for the fitted parameters, for four galaxies in the sample. Contours of stellar mass vs. age are shown 
on the left, and E{B — V) vs. age on the right. The dark and light blue regions represent 68% and 95% confidence intervals respectively, 
and the red X marks the adopted best fit. The top two objects are characteristic of the young, low stellar mass galaxies in the sample; the 
third row shows a typical galaxy, with fit very well-constrained by the addition of mid-IR IRAC data; and the bottom panel shows one of 
the most massive galaxies in the sample, which can only be fit by an old stellar population. Plots are shown with the same axes for ease 
of comparison. 



For 11 additional objects we measure an upper limit on 
the velocity dispersion because the measured line width 
is not significantly higher than the instrumental resolu- 
tion, as described below. 12 objects have FWHM less 
than the instrumental resolution due to noise, and 6 
have been rejected due to significant contamination from 
night sky lines. We calculate the velocity dispersion 
a = FWHM/2.35, where FWHM is the fuU width at 
half maximum after subtraction of the nominal instru- 
mental resolution (15 A in the K-band) in quadrature. 
Uncertainties (Acup and Acrdown) are determined by cal- 
culating (Tniax and cTijiin from the raw line width W and 
its statistical uncertainty AW; W + AW gives (Tmax and 
W — AW gives cTmin- This results in a larger lower bound 



on the error when the line width is close to the instrumen- 
tal resolution, and an upper limit on a when W — AW 
is less than the instrumental resolution. In such cases 
we give one standard deviation upper limits of cr -|- AtJup . 
As we discuss further below, we have also measured the 
spatial extent of the Ha emission for each object; for 14 
objects, the FWHMhq is less than the slit width. In such 
cases the effective resolution is increased by the ratio of 
the object size to the slit width; we have made this cor- 
rection in the calculation of a for these 14 objects. The 
correction is usually < 10%. 

Velocity dispersions are given in column 6 of Table 01 
From the 85 measurements we find {a) — 120 ±9 km s~^ 
(108 ± 5 km s~^ with AGN removed), with a standard 
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Fig. 5. — Histograms showing the distributions of the results of 
the SED modehng. From left to right and top to bottom, we show 
E{B — V), age, star formation rate, and stellar mass. Statistics of 
the distributions are given in the text. 




Fig. 6. — A comparison of the stellar masses of the full sample 
of 461 UV-selected galaxies (red) and the 93 galaxies for which 
we have model SEDs and Ha spectra (blue). The mean stellar 
mass of the full sample is 3.0 ± 0.5 X 10^'^ Mq, while that of the 
Ha subsample is 3.7 ± 0.5 X 10^'' Mq. For the purposes of this 
comparison we use constant star formation models for all objects, 
and the K magnitudes of the NIRSPEC sample have not been 
corrected for Ha emission. 

deviation of 86 km s~^ , while counting the upper hni- 



o-l 



(101 ± 5 



^ without AGN) with a standard deviation of 85 
^ . The observed hne widths could be caused by 



its as detections gives (a) = 112 ± 9 km s 
km s 
km s 

random motions, rotation, merging, or, most likely, some 
combination of all of these. Galactic outflow speeds (as 
measured by the velocity offsets between Ha, Lya, and 
the rest-frame UV interstellar absorption lines) show no 
correlation with the line widths, and the AGN fraction of 
~ 4% means that broadening by AGN is usually not sig- 
nificant. Dynamical masses can be calculated from the 
line widths via the relation 



The factor C depends on the galaxy's mass distribution 
and velocity field, and may range from C < 1 to C > 
5. The value of C depends on the mass density profile, 
the velocity anisotropy and relative contributions to a 
from random motions or rotation, the assumption of a 
spherical or disk-like system, and possible differences in 
distribution between the total mass and that of tracer 
particles used to measure it. Obviously the definition of 
the radius r is also crucial. Most of these factors are 
unknown for the current sample. 

Under the assumption that a disk geometry is appro- 
priate for gas-rich objects (see ^, we begin with the re- 
lationship Mdyn = Wtruo^/C- We incorporate an average 
inclination correction (wtruc) — 7''/2(wobs), and for Wobs 
use the observed velocity half-width (uobs = FWHM/2 = 
2.35(t/2). Combining the constants and writing the equa- 
tion in terms of a then gives C ~ 3.4. For the galaxy size 
fHa we use half of the spatial extent of the Ha emission, 
after deconvolution of the seeing. The galaxies are spa- 
tially resolved in almost all cases; for those few that are 
not, we use the smallest size measured under the same 
seeing conditions as an upper limit. We find a mean 
and standard deviation (rua) = 0.7 ± 0.3" {^ 6 kpc), 
approximately the same as the typical isophotal radius 
used for photometric measurements. Because it is diffi- 
cult to continuously monitor the seeing while observing 
with NIRSPEC, some uncertainty is introduced by our 
inexact knowledge of the seeing during each observation. 
This uncertainty is typically 0.1-0.2", but for those ob- 
jects that are near the resolution limit it may be much 
larger. Because of this issue, the sizes are uncertain by 
~ 30%. Note that most previous calculations of dynam- 
ical mass at high redshift have used C = 5 (appropriate 
for a uniform sphere; e.g. iPett ni ^1^ ^l]l2001t l Erb et all 



^£jPg£t 
il.120041: 



M, 



dyn 



Ca^r/G 



(1) 



[2003'. "2004*: 'Shaplev et al.l2004Hvan Dokkum et alJl2004t 
l^inbank et al. 200^^ 

Values of thq are given in column 5 of Table 21 and 
dynamical masses in column 8. Using the 85 galaxies 
with well-determined a, we find a mean dynamical mass 
(Mdyn) = 1-2 ± 0.4 X 10" Mq. The two largest dynam- 
ical masses in the sample are Mdyn = 5.7 x 10^^ Mq 
(Q1700-MD174) and Mdyn = 3.8 x lO^^ M© (Q1700- 
MD94); both of these objects are AGN, and their broad 
Ha lines are therefore not a reflection of the gravita- 
tional potential. Neglecting these two objects, and the 
other AGN which have Mdyn -- 6 - 9 x 10^° Mq (HDF- 
BMZ1156, Q1623-BX151 and Q1623-BX663), results in 
(Mdyn) = 6.9 ±0.6 X 10^° Mq, with a standard deviation 
of 5.8 X 10^° Mq. Approximately 25% of the non-AGN 
have dynamical masses Mdyn > 10^"'^ M©, with Q1623- 
BX376 (Mdyn = 3.7 x 10" M©) the largest. 

Uncertainties in the dynamical masses are probably 
dominated by the constant C, which may vary by a 
factor of a few. Our line widths may also occasionally 
suffer from contamination by night sky lines; we have 
removed the most obvious cases, but repeated obser- 
vations with varying signal-to-noise ratios (S/N) show 
that this can occasionally affect the line profile by up to 
^ 30%. Complex spatial and kinematic structure can 
also affect the observed vel ocity dispersion, as discussed 
previously llErb et alJl2003|) and recently emphasized by 
iColina et al.ll)200,'T() . The Ha emission also may not trace 
the full potential, especially in low S/N spectra of high 
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redshift objects. Even in some local starburst galaxies 
the Hg emission does not fu lly sample the rotation curve 
(jLehnert fc HeckmarJll99q) . Indeed, the mean dynam- 
ical mass we derive here is ~ 15 times lower than the 
typical halo mass as inferred from the galaxies' cluster- 
ing properties (Ad elbereer et al. 2005b'l . 

4.1. Comparisons with Stellar Masses 

Turning next to a comparison of the galaxies' dynam- 
ical and stellar masses, in Figure we plot dynamical 
mass vs. stellar mass for the 68 galaxies for which we 
have determined both. The dashed line shows equal 
masses, while the dotted lines on either side indicate 
approximate uncertainties of ±0.8 dex (a ~ 50% un- 
certainty in both masses, added in quadrature; this in- 
cludes systematic uncertainties in the dynamical mass 
model as well as the measurement uncertainties shown 
by the error bars in the lower right corner). Exclud- 
ing AGN, the mean dynamical mass is '^ 2 times higher 
than the mean stellar mass. Significant correlation be- 
tween stellar and dynamical m ass is obse rved in the local 
universe fe.g. iBrinchmann fc Ellis 2003), but such a cor- 
relation will exist only if the stellar mass makes up a 
relatively constant fraction of the dynamical mass over 
the full range of stellar masses. If low stellar mass galax- 
ies have large gas fractions , as is seen in local g alaxies 
l|McGaugh fc d e Blok 1997; Bel l fc de Jo ng 2000), they 
may not have correspondingly small dynamical masses. 
With these caveats in mind, we perform Spearman and 
Kendall r correlation tests (neglecting the AGN) and 
find a probability P = 0.01 that the masses are uncorre- 
lated, for a significance of the correlation of 2.5(T. This 
is consistent with our qualitative expectations from Fig- 
ure 13 where most of the points lie well within the dot- 
ted lines marking a factor of ^ 6 difference between the 
masses (though we see a significant number of outliers 
with Mdyn ^ M^,, to be discussed further below). 

Approximately 20% of the galaxies in the sample have 
Ali, > Afdyn, clearly an unphysical situation. The dis- 
crepancies are usually a factor of 2-3 and always less 
than a factor of ~ 6, which can be easily explained by 
the factor of a few uncertainty in both masses. In any 
case, because the Ha linewidths may not sample the full 
potential it would not be a surprise to find some objects 
with M^ > Mdyn- 

It is potentially more interesting to consider the small 
but significant fraction of the objects for which A/dyn ^ 
Mi,; we focus on the r^ 15% of the sample (11/68) for 
which Mdyn/M* > 10 (one of these is the AGN Q1700- 
MD94, which has a dynamical mass of 3.8 x 10^^ M© 
due to its very broad Ha line; this is clearly a prob- 
lematic estimate of dynamical mass, and we exclude this 
object in the following discussion). These objects are 
marked with open diamonds in Figure |7| and in many 
of the subsequent figures as well. The large differences 
between stellar and dynamical masses suggest the follow- 
ing possibilities: 1) these are galaxies that have recently 
begun forming stars, and thus have small stellar masses 
and large gas fractions; or 2) we significantly underes- 
timate the stellar mass for up to ~ 15% of the sample, 
presumably because of an undetected old stellar popula- 
tion. One should also consider the possibility that con- 
tamination from winds, shock ionization or AGN causes 
an overestimate of the dynamical masses of ^ 15% of 



the sample. There is no evidence for this third possi- 
bility, however; except in the few cases of known AGN, 
none of the objects show evidence of ionization from any 
source other than star formation, either from their UV 
spectra or from high [N II] /Ha ratios. The obj ects with 
Afdvn ^ Mi, have low stellar masses, and as lErb et alJ 
l)2006aj) show in the context of estimating metallicity as 
a function of stellar mass, [N II] is not detected even in 
a composite spectrum of 15 such objects. 

The best fit ages and Ha equivalent widths IjErb et alJ 
(|2006Ij) discuss the equivalent widths more fully ) of the 
set of galaxies with Mdyn/Af* > 10 favor the first pos- 
sibility. As shown in Figure |H1 there is a strong corre- 
lation between M^yn/Mi, and age, and, with somewhat 
more scatter, between Mdyn/Mi, and Whq (the corre- 
lations have 5 and 4 a significances respectively). Ob- 
jects with Mdyn ^ Mi, have young best fit ages, and 
tend to have high Ha equivalent widths as well, with 
200 < Wiia < 1300 A. Note that these are not correla- 
tions of independent quantities; age and stellar mass are 
correlated by definition in the constant star formation 
models we use for the majority of the sample, and the 
equivalent width depends on the Ks magnitude, which 
is also closely related to stellar mass. This, combined 
with the fact that we see less dispersion in dynamical 
mass than stellar mass, accounts for the correlations in 
Figure El 

It is important to assess whether or not the Mdyn ^ 
Mi, galaxies are really significantly younger than the rest 
of the sample, given the well-known degeneracies between 
age and extinction in SED modeling and the consider- 
able uncertainties in the age of a typical object. Confi- 
dence intervals for the a,ge-E{B — V) fits for two of the 
Mdyn > Mi, galaxies (Q1623-BX502 and Q2343-BX493) 
are shown in Figure 01 These are representative of this 
set of objects. Young ages (and correspondingly high val- 
ues oi E{B — V)) are clearly strongly favored, but there 
is a tail of acceptable solutions extending to higher ages. 
This set of objects is unique in favoring such young ages. 
To further test the significance of the young ages, we 
have divided the sample into four quartiles by Mdyn/M^ 
(with 16, 17, 17 and 17 galaxies in the quartiles), and per- 
formed K-S tests on the age distribution of each quartile, 
making use of the Monte Carlo simulations described in 
S]31 These simulations perturb the colors of each galaxy 
according to its photometric errors, and compute the 
best-fit model for the perturbed colors. We conducted 
10,000 trials for each galaxy, and generated lists of the 
best-fit ages for all trials in each quartile (for 160,000 
or 170,000 trial ages in each quartile). A two sample 
K-S test on all pairs of the four quartiles finds that the 
probability P ci that the ages in the quartile with the 
highest Mdyn/ Mi, are drawn from the same distribution 
as the remaining quartiles. This test suggests that even 
though the ages of individual galaxies may not be well 
constrained, the younger average age of the subsample 
with Mdyn ^ Mi, is highly significant. 

The best-fit ages from the SED fitting represent the 
age of the current star formation episode, however, not 
necessarily the galaxy as a whole, which brings us to pos- 
sibility 2) above: is there a significant underlying pop- 
ulation of old stars in these objects which causes us to 
underestimate the stellar mass? We address this question 
through the two-component models described in 
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Fig. 7. — Dynamical vs. stellar mass for the 68 galaxies with measurements of both quantities. Excluding AGN, we find (M*) = 
3.5 it 0.5 X 10^'' Mq, while the mean dynamical mass is {M^y^) = 7.1 ± 0.7 X 10^'' Mq. The dashed line indicates equal masses. The 
error bars in the lower right corner show typical uncertainties, which for dynamical masses include only measurement errors in rjja and a. 
Uncertainties in dynamical masses are dominated by the unknown factors in the mass model as discussed in the text; incorporating these 
uncertainties, we estimate approximate combined uncertainties of ~ 0.8 dex, as shown by the dotted lines. Galaxies with A/^yn/Af* > 10 
are marked with open diamonds, and we discuss possible reasons for their discrepancy in the text. 

which fit the SED with the superposition of a young 
burst and maximally old population to estimate the max- 
imum stellar mass. The maximal mass models, which fit 
all the flux redward of the ii'-band with an old burst, 
and then fit a young model to the UV residuals, result 
in stellar masses ~5-30 times larger for those galaxies 
with Mdyn/Afj, > 10. If we use the maximal mass mod- 
els for these objects rather than the single-component 
models in the comparison with dynamical mass, we find 
(Mdyn/M^) — 3, as compared to (Afdyn/A/^) — 54 us- 
ing the single component models. However, as described 
these models are poorer fits and result in im- 



plausibly high average SFRs. A further difficulty is that 
the mean age of the young component in these mod- 
els is 3.5 Myr, approximately 1/20 the typical dynami- 
cal timescale for these galaxies. The more general two- 
component modeling, in which the total mass comes from 



a varying linear combination of a maximally old burst 
and an episode of current star formation, does not favor 
such extreme models; for the M^yn ^ M^, galaxies in 
question, the use of these more general models increases 
the average stellar mass by a factor of three. 

More powerful constraints on the ages and masses of 
the Mdyn ^ M^, objects come from the spectra of the 
galaxies themse lves. Usi n g corn posite spectra of the 
current sample, lErb et aLl lj2006aj) show that there is a 
strong trend between metallicity and stellar mass in the 
z ~ 2 sample, such that the lowest stellar mass galax- 
ies also have low metallicities. These objects likely have 
large gas fractions, accounting for their relatively low 
metal content (we discuss the derivations of the gas frac- 
tions further in i©. These results strongly favor the 
hypothesis that the large dynamical masses of these ob- 
jects are due to a large gas mass rather than a significant 
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Fig. 8. — Left: The mass ratio AI^y^/Mi, vs. the best-fit age from SED modeling, and right, M^y^/Mi, vs. Ha equivalent width Wna- 
Galaxies with high Mdyn/M* have young ages, and also tend to have high equivalent widths. Typical uncertainties are shown by the error 
bars in the corner of each plot. Symbols are as in Figure 171 

population of old stars. From both considerations of dy- 
namical mass and baryonic mass (gas and stars, ^, it 
seems clear that objects with a best-fit stellar mass of 
^* ^ 10^ ^Q ^re not in fact significantly less massive 
than the rest of the sample. 



4.2. Velocity Dispersion and Optical Luminosity 

The correlation between the luminosity of elliptical 
galaxies and their velo city dispersion a (the Faber- 
Jackson relation, L oc g'^ jFaber &: Jackso n"1976^ is well- 
established in the local universe. Many proposals for the 
origin of this correlati on rely on self-regula ting feedback 
processes; for example . iMurrav et al.ll)200,^ propose that 
momentum-driven winds from supernovae and radiation 
pressure regulate the luminosity from star formation at 
L oc (7^. Other theor ies (e.g. Bcgclman & Nath 2005; 
iRobertson et al.l 120051 and references therein) use feed- 
back from accretion onto a supermassive black hole to 
produce a correlation between luminosity or stellar mass 
and velocity dispersion as a byproduct of the relation- 
ship between black hole mass an d velocity dispersion, 
Mrh oc a'^ jFerrarese fc Merrittil2000t iGebhardt et all 
l2000t iTremaine et al.ll2002D . In the local universe, these 
correlations exist in relatively homogeneous, quiescent 
galaxies; with the present data, we can test the rela- 
tionships between luminosity, stellar mass and velocity 
dispersion at a time in which the galaxies are building 
up much of their stellar mass. 

Previous efforts to find such relations at high redshift 
have generally suffered from small sample sizes and re- 
sulted in w eak or absent correlations with large aniounts 
of scatter f Pettini et al."2001l: Ivan Dokkum et all 120041: 
iSwinbank et al..:20041 . We revisit the relation between 
a and optical luminosity with a sample of 77 z '^ 2 UV- 
selected galaxies, augmented by 21 SCUBA galaxies pre- 
sented bylSwinbank et al. (2004), fou r of the Distant Re d 
Galaxies (DRGs) discussed bv lvan Dokkum e^ al. (2004), 
and 16 Lyman-break galaxies at z ^ 3 IjPettini et al. 
I2OOII give velocity dispersions for nine galaxies observed ''''°''*'''' wavelengths. 



in K and (usually) J bv lShanlev et al.|l20011 and we have 
subsequently observed seven more). Our total sample 
thus consists of 118 galaxies with a mean spectroscopic 
redshift (z) = 2.38 and a^ = 0.37. 

Figure ini shows the extinction-corrected rest- frame V- 
band luminosity plotted against the velocity dispersion a] 
the luminosity is determined by multiplying the best-fit 
SED with the redshifted V transi nission curve, a r id the 
extinction correction employs the iCalzetti et al.l l)2000f) 
extinction law and the best-fit values of E{B — V) from 
the SED modeling. We have used photometry from 
iShaplev et all l|200H) to re-model the SEDs of the z - 3 
LBGs, for consistency with the z ~ 2 sample. We 
also use the extinction-corrected magnitudes My from 
ISwinbank et al. ( 2004) , and u se the values of Ay deter- 
mined from SED modeling bv lvan Dokkum et al.l l|2004fl 
to correct their absolute V magnitudes. The solid line 
shows the local i-band^ Faber-Jackson relation derived 
from 9000 S loan Digital Sky Survey (SDSS) galaxies by 
iBernardi et al. (2003). Assuming no change in slope, the 
local relation is offset from the median of the data by 
1.35 dex in luminosity, as shown by the dotted line. The 
data are strongly correlated; the Spearman and Kendall 
T tests give a probability P = 2 x 10~^ that the data 
are uncorrelated, for a significance of 4.1cr. An attempt 
to fit both the slope and the zeropint of the correlation 
shows that the slope is not well-constrained, and depends 
strongly on the fitting procedure and the relative uncer- 
tainties assumed for L and a; however, it is not incon- 
sistent with ^ 4, the observed slope of the local Faber- 
Jackson relation. 

Though there is strong evidence for a correlation be- 
tween L and cr at z > 2, the scatter is considerable and 
too large to be accounted for by observational errors, a 
varies by a factor of ~ 3-4 at a given luminosity, and 
for most values of a we see at least an order of mag- 

^ The slope of the Faber-Jackson relation is the same within the 
uncertainties in all bands, while the zeropoint decreases slightly at 
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Fig. 9. — Velocity dispersion vs. extinction-corrected rest-frame V-band luminosity. Red circles show the current sample, blue xs are 
Lyman break galax ies at 2 ~ 3, gre en stars are the DRGs presented by van Dokkum ct al. (2004), and the purple squares represent the 
SCUBA galaxies of lijwinbank et alj |2004|, filled squares are those classified as star-forming galaxies, and open squares represent AGN). 
The solid line shows the local Faber-Jackson relation, and the dotted line shows the local relation shifted upward by 1.35 dex in order to 
pass through the median of the data. Typical uncertainties are shown by the error bars in the lower right corner. 



nitude variation in L. 
of the scatter may be 



Some insight into the sources 
obtained by accounting for the 



large variation in the mass-to-hght rat io M /L observed 
at z ~ 2. As shown by Shaolcv ct al] l)2005j) for a simi- 
lar sample of UV-selected galaxies, the rest-frame optical 
M/L varies by a factor of ~ 80, with the highest values of 
AI/L approaching the typical values seen in local galax- 
ies; we find very similar results for the current sample. 
In Figure [Till we account for this variation by decreasing 
the luminosity of each galaxy in the z ~ 2 sample by 
the difference between its value of M/L and the typical 
(M/L)o of an elliptical gala xy in the local univer s e. We 
use (Af/L)o = 3, as found bv lPadmanabhan et "all (|2004f) 
for elliptical galaxies in the SDSS; the average shift in lu- 
minosity is a factor of 10, and the shift ranges from a fac- 
tor of 2 to a factor of 120. The figure approximates what 
would be observed at z '^ if each galaxy were fixed at its 
current stellar mass and velocity dispersion and faded to 



match present-day M/L values. This prediction is par- 
ticuarly relevant for the current sample because the clus- 
tering properties of the z ~ 2 galaxies indicate that they 
are the likely prog enitors of elliptical galaxi es like those 
seen in the SDSS f^Adelherger enill2005H) . The solid 
line again shows the local Faber-Jackson relation, which 
is a reasonable fit to many of the points; however, the 
large remaining scatter and significant number of out- 
liers are consistent with the expectation (from the high 
ga s fraction s discussed in ^ and the high SFRs found 
by lErb et a l. 2006b) that most of the galaxies are still 
actively building up stellar mass. 

For the current data, the most relevant prediction 
of the rel ationship betw e en L and a at high redshift 
is that of iMurrav et al.l l)2005j) , who propose a limit- 
ing, Eddington-like luminosity for starbursts, Lm — 
(4/gc/G')a^, where fg refers to the gas fraction. In this 
model, momentum-driven winds powered by radiation 
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Fig. 10. — The V-band luminosities of the 2: ~ 2 sample shifted by 
the difference between the observed and local mass-to-light ratios 
as described in the text, plotted vs. velocity dispersion a. The solid 
line shows the local Faber-Jackson relation. Typical uncertainties 
are shown by the error bars in the lower right corner. 
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Fig. 11. — Bolometric luminosity i^{ia-|-uv) calculated from 
the extinction-corrected Ha star formation rates as described in 
the text vs. velocity dispersion, compared with the model of 
IMurrav e t "all 1200^ for gas fractions fg = 1.0, 0.1, and 0.01 (solid, 
dashed and dotted lines respectively). These three lines also in- 
dicate reasonable uncertainties in the model. The local Faber- 
Jackson relation is shown by the dot-dash line at right. Typical 
uncertainties are shown by the error bars in the lower right corner. 



pressure and supernovae would expel much of the gas 
from the galaxy at luminosities above Lm- At Lm, feed- 
back moderates the SFR, and luminosity does not in- 
crease sigificantly beyond Lm- Lm represents the bolo- 
metric luminosity from star formation, so the plot of Ly 
vs. a discussed above is not the most appropriate test 
of this prediction. Instead we estimate the bolometric 
luminosities J^(ir+uv ) from the extinction-corrected Ha 
SFRs determi ned by |Er b ct al. (2006b). We invert the 
prescription of lBell et al. (^005), who relate the star for- 
mation rate to the bolometric luminosity, to produce the 



following relation: 

i(iR+uv)/L0 = 1.0 X IQio SFR/Mo yr-^. (2) 

Here Lir represents the total IR luminosity, while Luv 
is the light from unobscured stars. Before us ing this 
expres sion we have converted our SF Rs to the iKroupal 
l|2(¥)l IMF used bv lBellet al!* ("200^. We find a mean 
(L(iR+uv)) = 4.6 X 10"'^^ L0, slightly higher than the aver- 
age bolometric luminosity (^(ir+uv)) — 2.3 x 10^^ L© es- 
timated for a different sample of UV-sclcctcd z ^ 2 galax- 
ies from 24 /im observations by Rcddv ct al. ( 2006) . The 
two determinations of i(iR+uv) can be compared directly 
for 11 g alaxies common to b oth samples, as shown in Fig- 
ure 5 of lReddv et alJl)2006D : the two measurements agree 
within the uncertainties, with {\ogL^"^) = 11.46 ±0.27 
and i^ogL^^^) = 11.59 ±0.20. 

Figure 1111 shows the inferred bolometric luminosity 
-^(iR+uv) plotted against velocity dispersion, along with 
lines indicating Lm for gas fractions fg — 1.0, 0.1, and 
0.01 (from left to right, solid, dashed and dotted lines). 
We also show the local Faber-Jackson relation. Although 
the three Lm lines indicate differences in the gas frac- 
tion, they also represent the uncertainty in Lm- They 
are not sufficient to discriminate between gas fractions; 
the objects with the highest inferred gas fractions (the 
.^dyn 3> Mi, galaxies marked with open diamonds) do 
not fall clos est to the f„ = 1.0 lin e. Within these uncer- 
tainties, the IMurrav et alJ l)2005|) model appears to be a 
reasonable description of the data; the points are strongly 
correlated, with 4. Oct significance and P ~ 7 x 10^^, and 
the Lj\/ lines bracket most of them. 

The observed or inferred luminosity is a complex su- 
perposition of light from current star formation and 
from formed stellar mass (and from dust emission, in 
the case of bolometric luminosities), and the star for- 
mation rate likely varies on a much shorter timescale 
than CT. A somewhat less complicated comparison may 
therefore be between stellar mass and velocity disper- 
sion, quantities that might be expected to be related 
through the Mbh — cr relation and the .Afbuieo — Mb 
relation (Magorrian et al. 1998; Marconi & Hunt 200 
as well as through the comparison of stellar and dynam- 
ical masses discussed above. We plot stellar mass A/^ 
vs. CT in Figure El this is similar to Figured although 
here we add the z ~ 3 LBGs and the DRGs. The cor- 
relation has a significance of 3.6ct, about the same as 
the CT — L correlation; the low stellar mass galaxies with 
large velocity dispersions discussed above are clear out- 
liers (marked with open diamonds; they are included in 
the above correlation test). The dotted line is the local 
relation between bulge stellar mass and velocity disper- 
sion, c onstructed by cq r nbin ing the M\,^\^a-MBH rela- 
ti on of [M arconi fc HuntI i2003.) and the Mbh-o' relation 
of iTremaine et al.l l|2002() : it provides a surprisingly ac- 
curate upper envelope in ct on the plot, as with the ex- 
ception of AGN and the Mdyn ^ Mi, objects discussed 
above (marked with open diamonds), virtually all of the 
points lie on or to the left of this line. 

At least some theoretical predictions indicate that the 
correlation between M^, and ct should be strong at high 
redshift. The results of one such study are shown by 
the open symb ols, from the numerical simulations of 
iRobertson et al. (,2005. see their Figure 3). These simula- 
tions test the evolution of the Mbh — c relation with red- 
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shift using simulations of merging galaxies, incorporating 
feedback from black hole growth such that ^ 0.5% of the 
accreted rest mass energy heats the gas. The Mbh — cr 
correlation then resu lts from the regulation of black hole 
growth by feedback i|Di Matteo et al.ll2005^ . The simu- 
lation results for z = 2 (open squares) and z = 3 (open 
circles) describe the upper envelope in observed a reason- 
ably well, but the data exhibit considerably more scatter. 
Given the parameters of the simulations, this is not a sur- 
prise. The velocity dispersions in the models (which are 
the stellar velocity dispersions, not the gas velocity dis- 
persions we use) are measured well after the merger and 
accompanying AGN feedback, when the star formatio n 
rate has dropped to nearly zero l)Di Matteo et al.ll2005|) . 
while the data represent galaxies at a variety of stages 
of an earlier starburst phase. The simulation sample is 
therefore considerably more homogeneous than the data 
set. 

A large number of galaxies in our data set have low 
velocity dispersions for their stellar mass relative to the 
local correlation (or high stellar mass for their velocity 
dispersion). In fact this is the opposite of the trend pre- 
dicted by the simulations, which indicate weak evolution 
with redshift in the Mbh — cr and Mi, — a relations, in the 
sense that velocity dispersions increase at a given stellar 
mass at higher redshifts (attributed to the steeper po- 
tential wells of high redshift galaxies). This difference 
may reflect the fact that we are probably underestimat- 
ing the true velocity dispersion in many cases, since we 
are only sensitive to the highest surface brightness star- 
forming regions, which may not sample the full poten- 
tial. It will be interesting to see if deeper spectra with 
larger telescopes result in a higher average velocity dis- 
persion for this sample, and an accompanying reduction 
of the scatter. This plot also provides a hint of an evo- 
lutionary sequence: if the Mdyn/Mi, > 10 objects are 
indeed young galaxies with high gas fractions, they may 
evolve upward on this plot, significantly increasing their 
stellar mass with relatively little change in velocity dis- 
persion until they fall on or to the left of the dotted 
line along with the rest of the sample. It is also pos- 
sible that the velocity dispersions of these objects will 
decrease over time, bringing them closer to the local re- 
lation, if they los e a large fraction o f their gas to outflows 
(as suggested bv lErb et al.ll2006a|) . From the observed 
correlations of a, L, and M^,, and their accompanying 
scatter, we conclude that the processes which create cor- 
relations between these properties in the local universe 
are underway at z ^ 2. 

4.3. A Comparison with z ^ S Lyman Break Galaxies 

Close examination of Figure El shows that the velocity 
dispersions of the z ~ 2 sample and the z ^ 3 LBGs have 
different distributions. We compare the two directly in 
FigureEl which shows histograms of the velocity disper- 
sions of the z r^ 2 galaxies and the LBGs, after excluding 
AGN. The mean and the error in the mean of the z ~ 2 
sample is (a) — 108±5 km s~^ , while for the z ~ 3 LBGs 
it is (a) = 84 ± 5 km s~^ ; the difference between the two 
means is of 5a significance. The z ~ 2 distribution is sig- 
nificantly broader; objects with low velocity dispersions 
are common at z '^ 2 as well as at z ~ 3, but the lower 
redshift sample contains a large number of objects with 
(T > 130 km s^^ . Only one such galaxy is present in the 
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Fig. 12. — Stellar mass vs. velocity dispersion. The dot- 
ted line shows the local relation between bulge stellar mass and 
velocity dispersion, constructed by combining the M^^^i^^-Mgu 
relation of 'Marconi & Hunt" ^2003) and the Mbh^'^ relation of 
ITfemainc ct al. 12002). Open squares and circles show the results 
of the simulations of Roberts on et al, ^2005), which test the evolu- 
tion of the Mbh — <^ relation with redshift. Typical uncertainties 
are shown by the error bars in the lower right corner. 

3 galaxies have an average dy- 



3.6 ± 0.6 X 10^° Mq, a factor of 



z ~ 3 sample. The z 
namical mass (M^yn) = 
two smaller than that of the z '^ 2 galaxies; the average 
stellar mass of the LBGs is also a factor of '^ 2 lo wer than 
that of the z ~ 2 galaxies f Shaplev et al.ll2001 fl. though 
this comparison is complicated by sample selection ef- 
fects. As described in detail below, these differences in 
velocity dispersion and dynamical mass appear to reflect 
real differences between the z ~ 2 and z ~ 3 samples. 




80 Galaxies 
<z) = 2.23±0.23 
<<7) = 108±41 km s 

30 Galaxies 
<z> = 3.11±0.14 
(a) = 84±28 km s" 



a (km s"') 

Fig. 13. — The distributions of the velocity dispersions of 2 ^ 2 
(larger blue histogram) and z ~ 3 (smaller red histogram) galaxies. 
At 2 ~ 2, (T is derived from the width of the Ha emission line, while 
[O III] is used for the 2^3 sample. 

Possible selection effects are a concern in attempting 
to understand the different distributions of a in the two 
samples. The z ~ 3 line widths are measured from 
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[O III]A5007, while Ha is used for the z ^ 2 galaxies. 
To see if systematic differences between line widths mea- 
sured from the two lines might be responsible for the 
effect, we have measured line widths from [O III]A5007 
for 8 of the z ^ 2 galaxies which also have Ha line widths. 
The eight galaxies have Ha velocity dispersions ranging 
from 60 to 150 km s~^ . There is no systematic differ- 
ence in the velocity dispersion from the two lines, with 
(chq) = 95 km s~^ , and (ctjoiii]} = 97 km s~^ . The 
individual values of a usually agree well, and when they 
do not, contamination by sky lines, particularly in the 
iJ-band where we measure [O III], appears to be the 
source of the discrepancy. There is therefore no evidence 
that the use of different lines is a signifcant issue in the 
comparison of the two samples. 

Another possible factor is the higher redshift of the 
LBGs. Surface brightness is a strong function of redshift, 
and it is possible that we are seeing only the bright- 
est, most central regions of the LBGs compared with 
flux from a larger region of the z ~ 2 sample. If this 
is true, we might expect to observe less emission from 
the faint, low surface brightness portions of the galax- 
ies at z ~ 3 than at z ~ 2; in other words, the z ~ 3 
galaxies may appear smaller, beyond the expected scal- 
ing in size with redshift fe.g. lFerguson et al.ll2004|) . We 
have measured the spatial extent of 18 of the [O III] 
emission lines, for comparison with the Ha sizes; as in 
the z ~ 2 sample, most of the lines are spatially re- 
solved, and (rjoiii]) = 0.6"~ 4.5 kpc. The mean value 
of rjoiii] is slightly smaller than the mean Ha spatial 
extent, for which we find (rHa) = 0.7"~ 6 kpc. This 
difference is almost identical to the mean size difference 
fou nd between UV - bright galaxies at z '^ 3 and z '^ 2.3 
by iFerguson et al.l l)2004D . who find that galaxy radii 
scale with redshift approximately as the Hubble parame- 
ter H~^{z), in accordance with theoretical expectations. 
In our adopted cosmology, this scaling predicts that the 
z ~ 3 LBGs should have mean radii ^ 70% smaller than 
the z ~ 2 sample, entirely consistent with our obser- 
vations. We therefore have reason to believe that the 
observed size differences reflect real differences between 
the two samples, rather than surface brightness effects. 

The selection criteria for observation with NIRSPEC 
were somewhat different for the two samples. The z ~ 
3 galaxies were chosen primarily because of their UV 
brightness, with additional objects selected because of 
their proximity to a QSO sightline. Unlike some of the 
observations of z '^ 2 galaxies, objects with bright or 
red near-IR magnitudes or colors were not favored. If 
such bright or red objects make up most of those with 
large line widths, this could be a plausible explanation 
for the differences. An examination of the list of z ^ 
2 galaxies with a > 120 km s~^ shows that some of 
them were indeed selected for their bright or red near- 
IR properties, but an approximately equal number were 
not. The velocity dispersions of z '^ 2 galaxies along 
QSO sightlines, which should be representative of the 
sample as a whole, have the same distribution as that 
of the full z ~ 2 sample, indicating that the selection of 
near-IR bright or red objects is probably not responsible 
for the presence of objects with large velocity dispersions 
in the z ^ 2 sample and their absence among the LBGs. 

Another possibility is that [O III] is more difficult to 



detect in massive galaxies; with increasing metallicity, 
the ratio of [O III]/[0 II] decreases, as does the ratio of 
the [O III] and [O II] lines to H/3 (or Ha) for metal- 
licities above ~ 1/3 .Zq. Thus the relative fluxes of 
Ha and [O III] might be expected to depe nd strongly 
on metallicity and hen ce stellar mass (jTremonti et alJ 
120041 lErb et al.ll2006a(l . Although uncertainties in flux 
calibration make the ratios uncertain, in our sample of 
z '^ 2 galaxies with both Ha and [O III] measurements 
FHa/FiQiii] ranges from ^^ 0.5 to ~ 5 as stellar mass 
increases from ~ 10^ to ^ 10^^ M©, and [O III] is not 
detected for one of the most massive galaxies (we do de- 
tect H/3 for this object). We see that [O III] is weak 
in massive, relatively metal-rich galaxies, and this could 
contribute to the absence of objects with large line widths 
in the z - 3 NIRSPEC samples^. Not aU of the objects 
with large velocity dispersions in the z ~ 2 sample have 
large stellar masses, however; ~ 60% of the objects with 
Afdyn ^ M^, havc a > 130 km s^^ , and these low stellar 
mass objects also have low metallicities which should re- 
sult in strong [O III] lines. The absence of such galaxies 
from the z '^ 3 sample (there may be one; we have no 
near-IR photometry and no stellar mass estimate for the 
z '--^ 3 galaxy with the largest velocity dispersion) is un- 
likely to be explained by selection effects, and probably 
indicates that low stellar mass galaxies with large line 
widths are rarer at z ~ 3 than at z ^ 2. 

There are reasons to expect physical differences be- 
tween the z '--^ 2 and z ~ 3 samples. We have already 
noted that galaxies at z '^ 3 are seen to be more com- 
pact than those at z '^ 2, in agreement with expectations 
from hierarchical galaxy formation theory. An analysis 
of the correlation lengths of the LBGs and the z ~ 2 
sample indicates that the z ~ 2 galaxies reside in ha- 
los ^ 3 times m ore massive than those hosting the LBGs 
IjAdelberger et al. 2005b) . The smaller dynamical masses 
of the z ~ 3 galaxies are in good agreement with this 
result. We conclude that the lower average velocity dis- 
persion of the LBGs is not fully explained by selection 
effects, and that the differences in a between the z ~ 2 
and z '^ 3 samples are likely to indicate real differences 
in the typical masses of galaxies in the two samples. 

5. SPATIALLY RESOLVED KINEMATICS 

In our initial study of 16 Ha spectra of z '^ 2 galax- 
ies, nearly 40% of the sample showed spatially resolved 
and tilted emission lines (Erb ct al. 2003). In the cur- 
rent, enlarged sample of 114 objects the fraction is much 
smaller, 14/114 or 12%; this is probably both because of 
exceptionally good conditions during our first observing 
run (we observed 6 of the 14 objects during this run) 
and simply because of small number statistics. Velocity 
shear may be caused by rotation, merging, or some com- 
bination of the two, and whether or not it is detected 
in a nebular emission line depends on the size, surface 
brightness and velocity structure of the object, its incli- 
nation, and the alignment of the slit with respect to the 
major axis. The seeing during th e obser vations plays a 
crucial role as well, as lErb et al] l)2004() show with re- 
peated observations of the same object. Inclinations and 

* In principle we would still recognize z ~ 3 objects with large 
velocity dispersions via their broad H/3 lines, but in many cases H/3 
is too weak for a reliable measurement 
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Fig. 14. — Observed velocity as a function of slit position for objects with spatially resolved and tilted Ha emission lines. The seeing for 
these observations was ~ 0.5", so the points shown are highly correlated, with approximately four points per resolution element. We plot 
one point per pixel to show the observed velocity field clearly, but the blurring of the seeing means that these diagrams do not represent 
the true velocity structure of the galaxies. 



major axes are unknown for nearly all of our sample, as 
is the primary cause of the velocity shear, and therefore 
this fraction of 12% should be considered a lower limit to 
the fraction of rotating or merging objects in the z ^ 2 
sample. The true fraction cannot be determined given 
the limitations imposed by the seeing. 

In order for the tilt of an emission line to be consid- 
ered significant, we required that the observed spatial 
extent be at least 1.5 times larger than the seeing disk, 
and that the peak-to-peak amplitude of the shear be at 
least four times larger than the typical velocity uncer- 
tainty. Five of the 14 objects with spatially resolved ve- 
locity shear have not been p revious ly discussed (six are 
descri bed bvlErb et alJl20na two bv lErb et al1l2n04 and 
one bvlShaol ev et al.|l2004D . We plot the observed veloci- 
ties with respect to systemic as a function of slit position 
for these five objects in FigurelT ^ we also include Q1 623- 
BX528, the object discussed bv lShaplev et aLll2nf)4i be- 
cause such a diagram was not presented in that paper. 
These diagrams are constructed by stepping along the 
slit pixel by pixel, summing the flux of each pixel and 
its neighbor on either side to increase the S/N, and mea- 
suring the centroid in velocity at each position. We em- 
phasize that these are not rotation curves of the variety 
plotted for local galaxies. The points are highly corre- 
lated because of the seeing (~0.5" for the observations 
presented here, thus there are '^ 4 points per resolution 
element; we plot one point per pixel to give a clear pic- 
ture of the observed velocity field), and recovery of the 
true velocity structure requires a model for the struc- 



ture of the obj ect because of th e degeneracies created by 
the seeing Csee lLaw et al.ll200g for a demonstration). Al- 
though the physical scale corresponding to the angular 
scale is given at the top of each panel for reference, the 
implied mapping between physical radius and velocity is 
not a true representation of the structure of the object. 

Although the seeing prevents us from determining the 
true velocity field of the objects, we can at least deter- 
mine a lower limit on the amplitude of the velocity shear 
from the uncorrelated endpoints of each curve. We cal- 
culate the observed Vc = (wmax — Vmin)/'^ for each ob- 
ject, where Umax and iJmin are with respect to the sys- 
temic redshift. This estimate is almost certainly less 
than the true velocity shear because the inclinations of 
the galaxies are unknown, in most cases we have made 
no attempt to align the slit with the major axis, and the 
seeing effectively reduces the velocities at the endpoints 
by mixing the light from the edges of the galaxy with 
emission from higher surface brightness, lower velocity 
regions toward the center. Even under ideal conditions, 
the Hq emission may not tr ace the true circular velocity; 
iLehnert fc HeckmanI l)199fi|^ find that in local starburst 
galaxies the regions of high Ha surface brightness sample 
only the inner, solid-body portion of the rotation curve. 
The observed values of Vc for the 14 objects with velocity 
shear are given in Table 01 and in the left panel of Fig- 
ure llSl thev are plotted against the velocity dispersion a. 
The dotted line marks equal values; we see that Vc ^ cr 
for most objects, though with considerable scatter. 

The relationship between a and the terminal velocity 
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Vc (we use Vc to refer to our observed velocity shear, and 
Vc to indicate the true circular velocity of a disk ) has been 
quantified in samples of more local galaxies. 'Ri x et aP 
|l997) find a ~ 0.614, while ^izzclla ct al. (200^ find 
Vc = 1.32 CT + 46 for elliptical and high surface brightness 
disk galaxies, for cr > 50 km s~^ and with velocities 
in km s~^ . It is clear that a and Vc of the galaxies in 
our sample do not follow these distributions; instead we 
find a mean {(j/vc) ^ 1.2 (while it is not possible for 
Vc to be less than cr, in the limit as the galaxy becomes 
spatially unresolved the observed Vc — > 0, and a/vc — > 
oo). Because our measurements of a are less affected 
by the seeing, they are undoubtably more reliable than 
our measurements of Vc- If we assume for the moment 
that our galaxies are rotating disks, we can use one of 
the local relations to predict Vc f or our sample. If th e 
z ~ 2 galaxies obey the relation of lPizzella et alJ l)2005j) . 
our measurements of Vc underestimate the true circular 
velocities by an average factor of ~ 2, and the full Ha 
sample has an average (Vc) ^ 190 km s~^ . We have 
already shown that a change in the seeing from ~ 0.5" 
to ~ 0.9" can reduce the observed fr. by a factor o f ~ 2 
while changing a by less than 10% l)Erb et al .120041) . so it 
is not unreasonable to assume that a change in resolution 
from ~ 0.1" or better to ~ 0.5" might have a similar 
effect. Deep observations of these objects at high angular 
resolution will be required to obtain a true measure of 

Vc. 

Next we ask whether or not the galaxies which display 
velocity shear are different from the rest of the galax- 
ies in any significant way. Figure [T31 shows that galaxies 
with shear span nearly the full range in velocity disper- 
sion; the mean value of a for the galaxies with shear 
is 119 km s~^ , while the mean a of galaxies without 
shear is 120 km s~^ . There is mild evidence that galax- 
ies with shear tend to be older; the median age of the 
galaxies with shear is 1434 Myr, compared to 509 Myr 
for the galaxies without shear. The galaxies with shear 
also have slightly higher stellar masses, with a median of 
4.4 X 10^° M© compared to 1.8 x 10^° for those without 
shear. A larger sample of galaxies with shear, and more 
uniform observing conditions, are required to determine 
whether or not these differences are significant. Assum- 
ing for the moment that they are, they may suggest that 
the rotation of mature, dynamically relaxed galaxies is a 
more important contribution to our observed shear than 
merging, which should not have a preference for older, 
more massive galaxies. 

Nothing in the data is inconsistent with ordered rota- 
tion as the primary cause of the observed shear. With 
this in mind we draw attention to the position-velocity 
diagram of Q2343-BX389, shown in the lower left panel of 
Figure ITU which appears to turn over or flatten at both 
ends in the manner of locally observed rotation curves. 
This is one of several objects from the 2 ~ 2 sample re- 
cently observed with SIN FONI, the near-IR integral fiel d 
spectrograph on the VLT IjForster Schreiber et al.''2006^ . 
The SINFONI position-velocity diagram also shows evi- 
dence for fiattening, but at higher velocities; the velocity 
difference is probably due to the ~ 50° misalignment of 
the NIRSPEC slit with the kinematic major axis and the 
deeper SINFONI integration which reveals flux at larger 
radii. Several but not all of the galaxies observed with 
SINFONI show signatures of rotating disks, and the role 
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Fig. 15. — A comparison of the velocity dispersion a and the 
observed velocity shear Vc = (iimax — i'min)/2 for the 14 objects 
with tilted lines. The dotted line shows equal values. We find 
(a/vc) = 1.2. 

of merging in z ~ 2 galaxies still remains unclear. Given 
its expected importance at high redshift it is likely that 
it makes some contribution to our observed velocities. 
Improved disentangling of these effects must await high 
angular resolution spectroscopy with the aid of adaptive 
optics; such instruments will be able to map the velocity 
fields of high redshift galaxies at a reso lution impossible 
with our current data (jLaw et al.ll200(T() . 

6. GAS AND BARYONIC MASSES 

We have used the stellar and dynamical masses deter- 
mined above to infer that the galaxies in our sample have 
significant gas fractions. It would be highly desirable to 
test this directly, for example with CO measurements, 
but such data would require extremely long integration 
times with current technology, making the assembly of 
a large sample of gas masses impossible. We can, how- 
ever, exploit the correlation between star formation and 
gas density to obtain an estimate of the gas masses. In 
star-forming galaxies in the local universe, the surface 
densit ies of sta r form ation and gas are observed to fol- 
low a ISchmidtl ljl959(l law, Ssfr = ^S^^, over more 
than six orders of magnitude in Ssfr l)Kennicuttl ll998^ . 
This empirical relation is usually explained by a model 
in which the SFR scales with density-dependent gravita- 
tional instabilities in the gas. The correlation has not yet 
been tested at high redshift because of the lack of mea- 
surements of gas masses, although the one well-studied 
example, the lensed z = 2.7 LBG MS1512-c B58, appears 
to be consistent with the local Schmidt law IjBaker et alJ 
I2004|) . Assuming that the galaxies in our sample obey 
such a law, we can use the SFRs and the galaxy sizes thq 
measured from the spatial extent of the Ha emission to 
compute their star formation densities, and thus esti- 
mate their gas densities and ma sses. We use the global 
Schmidt law of lKennicut'^ l|1998D : 



SsFR = 2.5 X 10" 



■^gas 



1 Mq pc- 



1.4 



Mq yr ^ kpc 



(3) 
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in combination with the conversion from Ha luminosity 
to SFR from the same paper, 



SFR (Mq yr-i) 



L(Ha) 



1.26 X 10*1 erg s"! 



(4) 



to create an IMF-independent relation between our ob- 
served Ha luminosity per unit area and the gas surface 
density 



Sgas = 1.6 X 10- 



Shq 



erg s ^ kpc 



0.71 



Mo pc-2. (5) 



The radii used bv iKennicuttI l)1998|) to compute sur- 
face densities approximately coincide with the edge of 
the galaxies' Ha-emitting disks; for our surface densi- 
ties, we take an area equal to the square of the FWHM 
of the (continuum-subtracted when necessary) Ha emis- 
sion (rHct in Tabled}. For L(Ha) we use the extinction- 
corrected Ha luminosity, incorp orating the factor of 2 
aperture correction discussed bv lErb et al.l l)2006bD . We 
then take Mg^s = ^^gasJ^Ha ^^ ^^ estimate of the cold gas 
mass associated with star formation, and calculate the 
gas fraction /i = Mgas/(Mgas + M-^)- The removal of the 
IMF dependence from the Schmidt law (which assumes 
a Salpeter IMF in the conversion from Ha luminosity to 
SFR) facilitates comparison between gas, dynamical and 
stellar masses, as we discuss below. It is likely to be dif- 
ficult to measure the objects' sizes with accuracy, given 
their often clumpy UV continuum morphologies as seen 
in high resolution images fe.g. lErb et al.l200.'Ul2fl04l) : the 
Ha morphology is likely to be complicated as well. As 
noted in 21 the object sizes have typical uncertainties 
of ~ 30%, and the corrected Ha fluxes of the objects 
may be uncertain by up to a factor of ^ 2; therefore the 
gas masses of individual objects are uncertain by a fac- 
tor of ^ 3, and the typical fractional uncertainty in /x is 
'^50%. Systematic uncertainties due to the scatter in the 
Schmidt law itself (0.3 dex) are an additional source of 
error. We focus here on overall trends, which depend on 
large numbers of objects and so are better determined 
than the parameters of individual galaxies; for example, 
the error in the mean gas fraction of a subsample of ~ 15 
objects is 15% or less. 

The mean inferred gas mass is (Afgas) = 2.1 ±0.1 x 10^" 
Mq, slightly lower than the mean stellar mass of (M*) = 
3.6 ±0.4 X 10^" Mq. The distributions of the two masses 
are shown in the histograms in the upper left panel of 
Figure El where the narrower blue histogram shows the 
gas masses and the broader red histogram shows the stel- 
lar masses. The range of inferred gas masses is clearly 
smaller than the range of stellar masses; this is because 
the dispersions in star formation rate and size are smaller 
than the more than two orders of magnitude variation we 
see in stellar mass. The two masses are plotted against 
each other in the upper right panel of Figure 1161 where 
it is apparent that the Afdyn/Af* > 10 objects (open di- 
amonds) have large inferred gas masses, and are excep- 
tions to a general trend of increasing gas mass with in- 
creasing stellar mass. The absence of points in the lower 
left corner of this plot is probably a selection effect, as low 
mass galaxies with low star formation rates are unlikely 
to be detected in our K-h&nd images (they may fall below 
our TZ < 25.5 magnitude limit for selection as well), and 
we are also less likely to detect Ha emission in X-faint 



objects. We next plot the gas fraction ^ vs. stellar mass 
in the lower left panel of Figure El and /i vs. age in the 
lower right. The trends of decreasing gas fraction with 
increasing stellar mass and age are strong, supporting 
our hypothesis that the Afdyn ^ A^T* objects are young 
starbursts with high gas fractions (low mass galaxies with 
low gas fractions may exist, but are probably too faint to 
be detected by our survey). Local galaxies with low stel- 
lar masses are also observed to have higher gas fr actions 
lIMcGaugh fc de BlokTgg?: IBell fc de ,lol m^m]] . 

Other studies of galaxies at high and low redshift may 
provid e additional insight into these results. iReddv et alJ 
l)2006D have recently used 24 /im data from the Spitzer 
Space Telescope to infer the bolometric luminosities of 
both optically-selected galaxies similar to those consid- 
ered here and near -IR selected objects (Daddi et al. 200.4t 
iFranx et alJl2003|) . They assess the relationship between 
gas fraction and stellar mass, finding a similar strong 
trend which extends to massive, nearly passively evolv- 
ing galaxies with low gas fractions that are not selected 
by our UV criteria. They also consider the dust-to-gas 
ratio of local and high redshift galaxies, finding that at 
a given bolometric luminosity galaxies at z ~ 2 are ~ 10 
times less obscured by dust than local galaxies. Thus 
we infer that FIR emission is likely to be weaker in high 
redshift gas-rich galaxies than in their low redshift coun- 
terparts with similar bolometric luminosities. Galaxies 
that are detected at 850 ^m have so far provided the 
best candidat es for direct detection of gas in high red- 
shift galaxies: iGreve et al.l (|2P05) use CO measurements 
to estimate the gas masses of five submillimeter galaxies 
at z = 1-3.5, finding a median molecular gas mass of 
3.0 X 10^° Mq within a typical radius of ~ 2 kpc. This 
is only slightly higher than the typical gas masses we in- 
fer from the Schmidt law, implying that the largest gas 
masses we infer may be testable with current technology. 
Though the gas masses of the submm galaxies appear 
to be similar to those we find here, the factor of ~ 3 
smaller sizes of the submm galaxies imply a gas surface 
density an order of magnitude higher than the average 
value of our sample. Such comparisons must be viewed 
with caution, of course, given the very different methods 
of inferring the gas content and sizes of the two sam- 
ples; it is possible that a more direct measurement of the 
molecular gas distribution of the current sample would 
reveal smaller emitting regions and higher gas surface 
densities. 

We can now compare the estimated baryonic masses 
A^Tbar = Afgas + Af^, with the dynamical masses in Fig- 
ure El to see if the addition of the inferred gas mass 
improves the agreement between stellar and dynami- 
cal mass. We use the same axes as Figure [7| to fa- 
cilitate comparison, and again mark the line of equal 
masses (dashed line) and a factor of ^ 6 difference be- 
tween the two masses (dotted lines). The agreement 
between the masses is greatly improved with the addi- 
tion of the gas; there is a significant correlation (4tT, 
with probability P = 8 x 10^^ that the data are un- 
correlated), and the masses are within a factor of 3 for 
85% of the sample (objects with limits on a are not in- 
cluded, though most of these are consistent with the rest 
of the data, and those that are inconsistent are highly 
uncertain). After excluding AGN (open circles), we find 
(Afgas + A'U) = 5.8 ± 0.6 X 10^" Mq; the mean dynamical 
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Fig. 16. — The relation of the gas mass inferred from the Schmidt law to the stellar population parameters. Upper left: Histograms of 
the gas (blue) and stellar (red) masses, showing that the range in stellar mass is significantly broader than the range in gas mass. Upper 
right: Gas mass vs. stellar mass. Gas masses are relatively constant across the sample, but when the galaxies with M^y^/Mi, > 10 are not 
considered we see an increase in gas mass with stellar mass. Lower panels: The gas fraction strongly decreases with increasing stellar mass 
(left) and age (right). The solid line and accompanying shaded region in the lower left panel show the minimum detectable gas fraction set 
by our Ha flux detection limit, for the mean (solid line) and range (shaded region) of sizes observed. Symbols are as in Figure|7] and the 
error bars in the corner of each plot show typical uncertainties. 



mass (Mdyn) = 7.1 ± 0.7 x 10^" Mq is 1.2 times higher. 
There are only two objects for which M(jyn/Afbar > 10; 
these are the AGN Q1700-MD94 (upper right), and the 
galaxy Q1623-BX376, which has an anomalously large 
velocity dispersion that may be influenced by its compli- 
cated spatial structure ( Erb et al. 2003) . 

Although the estimates of gas, stellar and dynamical 
masses considered here carry significant uncertainties, 
we emphasize that they are obtained from independent 
quantities (Ha luminosity and size, multi-wavelength 
photometry, and line width respectively). Together they 
provide a remarkably coherent scenario in which rela- 

tivel y massive {M^„,^ -t- M^ ^■} ^^'^ ^©' ^'^'^ ^'1^1° '" ^'^^^ 
Mr:^ : lAdelberger et al.ll2005b|) star-forming galaxies have 
high gas fractions which decrease with age. The observed 
variation of metallicity across the sa mple also support s 
this model, as we describe separately l|Erb et alJl2006a|) . 



In this scenario, the stellar mass of a galaxy depends 
strongly on its evolutionary state, as well as on its total 
mass. The lowest stellar mass galaxies in the sample have 
low stellar masses because they are young and have con- 
verted or expelled only a small fraction of their gas; such 
objects have stellar masses > 100 times lower than the 
most massive galaxies in the sample, while their inferred 
baryonic and dynamical masses are only < 30 times lower 
than those of the most massive galaxies. 

7. CONCLUSIONS AND DISCUSSION 

The goals of this paper have been to examine the kine- 
matic properties of star-forming galaxies at z ~ 2 as re- 
vealed by their Ha spectra; to compare their dynamical 
masses with stellar masses determined from population 
synthesis modeling; to look for trends between kinematic 
and stellar population properties; and to estimate the 
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stellar and gas content of the galaxies. Our main conclu- 
sions are as follows: 

1. The sample has a mean Ha velocity disperson 
(a) = 101 ± 5 km s~^ , excluding AGN. From a and the 
spatial extent of the Ha emission we estimate dynami- 
cal masses, finding (Mdyn) — 7.1 ± 0.7 x 10^" M©, again 
excluding AGN. Stellar and dynamical masses agree to 
within a factor of '^ 6 for most objects, consistent with 
observational and systematic uncertainties. The masses 
are correlated with 2.5a significance. However, 15% of 
the galaxies have Mdyn > 10 x M^; these objects have 
low stellar masses, young ages and tend to have high Ha 
equivalent widths, suggesting that they are young galax- 
ies with large gas fractions. 

2. We combine our Ha results with LBGs at z ~ 3 and 
other samples from the literature, and find that rest- 



frame optical luminosity (corrected for extinction) and 
velocity dispersion are correlated with 4ct significance. 
The large scatter prevents a robust determination of the 
slope, but it is not inconsistent with the local relation 
L oc cr^. The high redshift galaxies have much higher 
luminosities at a given velocity dispersion than would 
be predicted by the local Faber-Jackson relation. The 
observed correlation and its accompanying large scatter 
suggest that the processes which produce the L-a corre- 
lation in the local universe are underway at z ~ 2-3. 

3. A comparison of the Ha sample with a similarly 
rest-frame UV-selected sample of Lyman break galaxies 
at z '^ 3 shows that the average velocity dispersion is 
~ 30% larger at z ~ 2 than at z '^ 3, and the average 
dynamical mass of the z '^ 2 galaxies is ~ 2 times larger 
than that of the z '-^ 3 LBGs. This difference proba- 
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bly reflects real differences in the average masses of the 
galaxies at z ~ 2 and z ~ 3. 

4. Fourteen of the 114 galaxies with Ha spectra, or 
12%, have spatially resolved and tilted emission lines. 
On average, the observed amplitude of the velocity shear 
Vc is approximately equal to the velocity dispersion a. 
If the galaxies are rotating disks and follow the local re- 
lations between the true circular velocity Vc and a, our 
observations underestimate Vc by an average factor of 
~ 2, and the galaxies have (Vc) ^ 190 km s^^ . However, 
other factors including galaxy-galaxy mergers, slit posi- 
tion angle, and the seeing FWHM are significant sources 
of uncertainty in Vc- 

5. Using the empirical correlation between star forma- 
tion rate per unit area and gas surface density, we esti- 
mate gas masses and gas fractions for the galaxies in the 
sample. The mean gas fraction is ~ 50%, and the gas 
fraction strongly decreases with increasing stellar mass 
and age. The combined mass in gas and stars is signifi- 
cantly better correlated with the dynamical mass than is 
the stellar mass alone. The sample spans a much wider 
range in stellar mass than dynamical mass (or inferred 
baryonic mass), indicating that there is considerable dis- 
persion in the fraction of gas that has been converted to 
stars by z '^ 2. 

These differing ranges in stellar and baryonic mass im- 
ply that galaxies of similar total masses began forming 
stars at different times, and that at a given redshift 
galaxies with a wide variety of total masses have just 
begun forming stars. For ~ 70% of the galaxies in the 
Ha sample, the best-fit ages imply formation redshifts 
2 < Zform ^ 4; most of these galaxies have baryonic 
masses between ^ 10^° and ~ 10^^ Mq. The remaining 
~ 30% have Zform > 4 (including the 15% of the sample 
which have best-fit ages equal to the age of the universe); 
these galaxies have IQi'^'' < Afbar < 10"* Mq. Mbar 
and Zform are strongly correlated across the full sample, 
with 5.4cr significance, while those with 2 < Zform ^ 4 
are correlated with 2.2a significance and there is no cor- 
relation between Mbar and Zform within the Zform > 4 
group. The two groups are shown in Figure 1171 where 
galaxies with Zform > 4 are shown by red points and 
galaxies with Zform < 4 are plotted in orange. Galax- 
ies with Mbar ;$ 10^^ Mq fall in both groups. These 
results are not unexpected, since total mass is not the 
only factor which determines when a galaxy begins to 
form stars; large scale environment and interactions with 
other galaxies are also significant. Galaxies in dense en- 
vironments will cross the critical density threshold for 



collapse earlier, and should therefore begin forming stars 
at earlier times. One of the fields in the present survey 
cont ains a signifi cant overdensity of galaxies at z = 2.3, 
and iSteidel et a l. (2005) show that galaxies within the 
overdensity are significantly older than average. It is also 
expected that mergers will trigger significant star forma- 
tion events; the mer ger-based numerical simulations of 
iSpringel et al.l l)2005l and references therein) are success- 
ful in reproducing many properties of galaxies. We can- 
not yet identify the likely triggers for star formation for 
galaxies on an individual basis, however; detailed mor- 
phological analysis and an improved quantification of the 
galaxies' environments will be useful in this regard. 

Because of uncertainties in the sizes and velocity dis- 
tributions of the z ~ 2 galaxies, the dynamical masses 
presented here have much room for improvement. These 
improvements will require deep observations of a large 
sample of objects, and in many cases high angular resolu- 
tion will be needed in order to discriminate between plau- 
sible models of the velocity field. Near-IR multi-object 
and integral field spectrographs, especially with adaptive 
optics, will make such observations feasible, and we an- 
ticipate that the galaxies presented here will be a useful 
sample from which to select targets for such observations. 
These kinematic measurements will provide otherwise 
unobtainable insights about the prevalence of mergers at 
high redshift and the growth of massive disk galaxies. We 
also anticipate studies at longer wavelengths which will 
improve our understanding of the galaxies' gas masses. 
All of these results will lead to a more refined picture of 
the assembly of mass at high redshift. 
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TABLE 1 
Galaxies Observed 



Object 


R.A. 


Dec. 


2Ha 


-R.^ 


G-n^ 


U„-G^ 


Ks^ 


J~Ks^ 


Exposure time (s) 


GDFb-BN88'^ 


00:53:52.87 


12:23:51.25 


2.2615 


23.14 


0.29 


0.68 






12 


X 720 


HDF-BX1055 


12:35:59.59 


62:13:07.50 


2.4899 


24.09 


0.24 


0.81 






2 


X 900 


HDF-BX1084 


12:36:13.57 


62:12:21.48 


2.4403 


23.24 


0.26 


0.72 






5 


X 900 


HDF-BX1085 


12:36:13.33 


62:12:16.31 


2.2407 


24.50 


0.33 


0.87 






5 


X 900 


HDF-BX1086 


12:36:13.41 


62:12:18.84 


2.4435 


24.64 


0.41 


1.09 






5 


X 900 


HDF-BX1277 


12:37:18.59 


62:09:55.54 


2.2713 


23.87 


0.14 


0.61 


21.26 


0.73 


3 


X 900 


HDF-BX1303 


12:37:11.20 


62:11:18.67 


2.3003 


24.72 


0.11 


0.81 


21.03 


2.31 


2 


X 900 


HDF-BX1311 


12:36:30.54 


62:16:26.12 


2.4843 


23.29 


0.21 


0.81 


20.48 


1.56 


4 


X 900 


HDF-BX1322 


12:37:06.54 


62:12:24.94 


2.4443 


23.72 


0.31 


0.57 


20.95 


2.16 


6 


X 900 


HDF-BX1332 


12:37:17.13 


62:11:39.95 


2.2136 


23.64 


0.32 


0.92 


20.68 


1.77 


3 


X 900 


HDF-BX1368 


12:36:48.24 


62:15:56.24 


2.4407 


23.79 


0.30 


0.96 


20.63 


1.81 


4 


X 900 


HDF-BX1376 


12:36:52.96 


62:15:45.55 


2.4294 


24.48 


0.01 


0.70 


22.13 


1.02 


4 


X 900 


HDF-BX1388 


12:36:44.84 


62:17:15.84 


2.0317 


24.55 


0.27 


0.99 


19.95 


1.78 


2 


X 900 


HDF-BX1397 


12:37:04.12 


62:15:09.84 


2.1328 


24.12 


0.14 


0.76 


20.87 


1.08 


3 


X 900 


HDF-BX1409 


12:36:47.41 


62:17:28.70 


2.2452 


24.66 


0.49 


1.17 


20.07 


2.31 


2 


X 900 


HDF-BX1439 


12:36:53.66 


62:17:24.27 


2.1865 


23.90 


0.26 


0.79 


19.72 


2.16 


4 


X 900 


HDF-BX1479 


12:37:15.42 


62:16:03.88 


2.3745 


24.39 


0.16 


0.79 


21.30 


1.57 


5 


X 900 


HDF-BX1564 


12:37:23.47 


62:17:20.02 


2.2225 


23.28 


0.27 


1.01 


19.62 


1.69 


2 


X 900 


HDF-BX1567 


12:37:23.17 


62:17:23.89 


2.2256 


23.50 


0.18 


1.05 


20.18 


1.23 


2 


X 900 


HDF-BX305 


12:36:37.13 


62:16:28.36 


2.4839 


24.28 


0.79 


1.30 


20.14 


2.59 


4 


X 900 


HDF-BMZ1156 


12:37:04.34 


62:14:46.28 


2.2151 


24.62 


-0.01 


-0.21 


20.33 


1.71 


4 


X 900 


Q0201-B13<= 


02:03:49.25 


11:36:10.58 


2.1663 


23.34 


0.02 


0.69 






16 


X 720 


Q1307-BM1163 


13:08:18.04 


29:23:19.34 


1.4105 


21.66 


0.20 


0.35 






2 


X 900 


Q1623-BX151 


16:25:29.61 


26:53:45.01 


2.4393 


24.60 


0.14 


0.88 






2 


X 900 


Q1623-BX214 


16:25:33.67 


26:53:53.52 


2.4700 


24.06 


0.39 


1.12 






4 


X 900 


Q1623-BX215 


16:25:33.80 


26:53:50.66 


2.1814 


24.45 


0.26 


0.53 






4 


X 900 


Q1623-BX252 


16:25:36.96 


26:45:54.86 


2.3367 


25.06 


0.07 


0.55 






3 


X 900 


Q1623-BX274 


16:25:38.20 


26:45:57.14 


2.4100 


23.23 


0.25 


0.89 






3 


X 900 


Q1623-BX344 


16:25:43.93 


26:43:41.98 


2.4224 


24.42 


0.39 


1.25 






2 


X 900 


Q1623-BX366 


16:25:45.09 


26:43:46.95 


2.4204 


23.84 


0.41 


1.03 






2 


X 900 


Q1623-BX376 


16:25:45.59 


26:46:49.26 


2.4085 


23.31 


0.24 


0.75 


20.84 


1.61 


4 


X 900 


Q1623-BX428 


16:25:48.41 


26:47:40.20 


2.0538 


23.95 


0.13 


1.06 


20.72 


1.02 


4 


X 900 


Q1623-BX429 


16:25:48.65 


26:45:14.47 


2.0160 


23.63 


0.12 


0.56 


20.94 


1.46 


2 


X 900 


Q1623-BX432 


16:25:48.73 


26:46:47.28 


2.1817 


24.58 


0.10 


0.53 


21.48 


1.76 


4 


X 900 


Q1623-BX447 


16:25:50.37 


26:47:14.28 


2.1481 


24.48 


0.17 


1.14 


20.55 


1.66 


3 


X 900 


Q1623-BX449 


16:25:50.53 


26:46:59.97 


2.4188 


24.86 


0.20 


0.61 


21.35 


1.88 


3 


X 900 


Q1623-BX452 


16:25:51.00 


26:44:20.00 


2.0595 


24.73 


0.20 


0.90 


20.56 


2.05 


3 


X 900 


Q1623-BX453 


16:25:50.84 


26:49:31.40 


2.1816 


23.38 


0.48 


0.99 


19.76 


1.65 


3 


X 900 


Q1623-BX455 


16:25:51.66 


26:46:54.88 


2.4074 


24.80 


0.35 


0.87 


21.56 


1.83 


2 


X 900 


Q1623-BX458 


16:25:51.58 


26:46:21.39 


2.4194 


23.41 


0.28 


0.85 


20.52 


1.23 


4 


X 900 


Q1623-BX472 


16:25:52.87 


26:46:39.63 


2.1142 


24.58 


0.16 


0.91 


20.80 


1.91 


4 


X 900 


Q1623-BX502 


16:25:54.38 


26:44:09.25 


2.1558 


24.35 


0.22 


0.50 


22.04 


0.99 


3 


X 900 


Q1623-BX511 


16:25:56.11 


26:44:44.57 


2.2421 


25.37 


0.42 


1.05 


21.78 




4 


X 900 


Q1623-BX513 


16:25:55.86 


26:46:50.30 


2.2473 


23.25 


0.26 


0.68 


20.21 


1.83 


2 


X 900 


Q1623-BX516 


16:25:56.27 


26:44:08.19 


2.4236 


23.94 


0.30 


0.82 


20.41 


2.29 


3 


X 900 


Q1623-BX522 


16:25:55.76 


26:44:53.28 


2.4757 


24.50 


0.31 


1.19 


20.75 


2.03 


4 


X 900 


Q1623-BX528 


16:25:56.44 


26:50:15.44 


2.2682 


23.56 


0.25 


0.71 


19.75 


1.79 


4 


X 900 


Q1623-BX543 


16:25:57.70 


26:50:08.59 


2.5211 


23.11 


0.44 


0.96 


20.54 


1.31 


4 


X 900 


Q1623-BX586 


16:26:01.52 


26:45:41.58 


2.1045 


24.58 


0.32 


0.87 


20.84 


1.79 


4 


X 900 


Q1623-BX599 


16:26:02.54 


26:45:31.90 


2.3304 


23.44 


0.22 


0.80 


19.93 


2.09 


4 


X 900 


Q1623-BX663 


16:26:04.58 


26:48:00.20 


2.4333 


24.14 


0.24 


1.02 


19.92 


2.59 


3 


X 900 


Q1623-MD107 


16:25:53.87 


26:45:15.46 


2.5373 


25.35 


0.12 


1.43 


22.43 


1.72 


4 


X 900 


Q1623-MD66 


16:25:40.39 


26:50:08.88 


2.1075 


23.95 


0.37 


1.40 


20.15 


1.74 


3 


X 900 


Q1700-BX490 


17:01:14.83 


64:09:51.69 


2.3960 


22.88 


0.36 


0.92 


19.99 


1.54 


3 


X 900 


Q1700-BX505 


17:00:48.22 


64:10:05.86 


2.3089 


25.17 


0.45 


1.28 


20.85 


2.17 


4 


X 900 


Q1700-BX523 


17:00:41.71 


64:10:14.88 


2.4756 


24.51 


0.46 


1.28 


20.93 


1.86 


4 


X 900 


Q1700-BX530 


17:00:36.86 


64:10:17.38 


1.9429 


23.05 


0.21 


0.69 


19.92 


1.23 


4 


X 900 


Q1700-BX536 


17:01:08.94 


64:10:24.95 


1.9780 


23.00 


0.21 


0.79 


19.71 


1.31 


4 


X 900 


Q1700-BX561 


17:01:04.18 


64:10:43.83 


2.4332 


24.65 


0.19 


1.04 


19.87 


2.43 


2 


X 900 


Q1700-BX581 


17:01:02.73 


64:10:51.30 


2.4022 


23.87 


0.28 


0.62 


20.79 


1.94 


2 


X 900 


Q1700-BX681 


17:01:33.76 


64:12:04.28 


1.7396 


22.04 


0.19 


0.40 


19.18 


1.27 


4 


X 900 


Q1700-BX691 


17:01:06.00 


64:12:10.27 


2.1895 


25.33 


0.22 


0.66 


20.68 


1.82 


4 


X 900 


Q1700-BX717 


17:00:56.99 


64:12:23.76 


2.4353 


24.78 


0.20 


0.61 


21.89 


1.49 


4 


X 900 


Q1700-BX759 


17:00:59.55 


64:12:55.45 


2.4213 


24.43 


0.36 


1.29 


21.23 


1.35 


2 


X 900 


Q1700-BX794 


17:00:47.30 


64:13:18.70 


2.2473 


23.60 


0.35 


0.58 


20.53 


1.37 


3 


X 900 


Q1700-BX917 


17:01:16.11 


64:14:19.80 


2.3069 


24.43 


0.28 


0.95 


20.03 


1.95 


3 


x900 


Q1700-MD103 


17:01:00.21 


64:11:55.58 


2.3148 


24.23 


0.46 


1.49 


19.94 


1.96 


900 -1- 600 


Q1700-MD109 


17:01:04.48 


64:12:09.29 


2.2942 


25.46 


0.26 


1.44 


21.77 


1.75 


4 


x900 


Q1700-MD154 


17:01:38.39 


64:14:57.37 


2.6291 


23.23 


0.73 


1.91 


19.68 


2.08 


3 


X 900 
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Object 


R.A. 


Dec. 


Zna 


■R'^ 


G-n'^ 


U^-G'^ 


Ks^ 


J -Ks^ 


Exposure time (s) 


Q1700-MD174 


17:00:54.54 


64:16:24.76 


2.3423 


24.56 


0.32 


1.50 


19.90 




4 


X 900 


Q1700-MD69 


17:00:47.62 


64:09:44.78 


2.2883 


24.85 


0.37 


1.50 


20.05 


2.60 


4 


X 900 


Q1700-MD94 


17:00:42.02 


64:11:24.22 


2.3362 


24.72 


0.94 


2.06 


19.65 


2.46 


3 


X 900 


Q2343-BM133 


23:46:16.18 


12:48:09.31 


1.4774 


22.59 


0.00 


0.19 


20.50 


0.66 


3 


X 900 


Q2343-BM181 


23:46:27.03 


12:49:19.65 


1.4951 


24.77 


0.12 


0.29 






4 


X 900 


Q2343-BX163 


23:46:04.78 


12:45:37.78 


2.1213 


24.07 


-0.01 


0.71 


21.38 


0.97 


4 


X 900 


Q2343-BX169 


23:46:05.03 


12:45:40.77 


2.2094 


23.11 


0.19 


0.72 


20.75 


1.05 


4 


X 900 


Q2343-BX182 


23:46:18.04 


12:45:51.11 


2.2879 


23.74 


0.14 


0.56 


21.60 


0.92 


4 


X 900 


Q2343-BX236 


23:46:18.71 


12:46:15.97 


2.4348 


24.28 


0.14 


0.71 


21.25 


1.60 


3 


X 900 


Q2343-BX336 


23:46:29.53 


12:47:04.76 


2.5439 


23.91 


0.40 


1.15 


20.80 


1.75 


4 


X 900 


Q2343-BX341 


23:46:23.24 


12:47:07.97 


2.5749 


24.21 


0.38 


0.89 


21.40 


2.07 


3 


X 900 


Q2343-BX378 


23:46:33.90 


12:47:26.20 


2.0441 


24.80 


0.26 


0.55 


21.90 


1.30 


4 


X 900 


Q2343-BX389 


23:46:28.90 


12:47:33.55 


2.1716 


24.85 


0.28 


1.26 


20.18 


2.74 


3 


X 900 


Q2343-BX390 


23:46:24.72 


12:47:33.80 


2.2313 


24.36 


0.24 


0.79 


21.29 


1.73 


4 


X 900 


Q2343-BX391 


23:46:28.07 


12:47:31.82 


2.1740 


24.51 


0.25 


1.01 


21.93 


1.35 


3 


X 900 


Q2343-BX418 


23:46:18.57 


12:47:47.38 


2.3052 


23.99 


-0.05 


0.37 


21.88 


1.76 


5 


X 900 


Q2343-BX429 


23:46:25.25 


12:47:51.20 


2.1751 


25.12 


0.30 


0.85 


21.88 


1.97 


4 


X 900 


Q2343-BX435 


23:46:26.36 


12:47:55.06 


2.1119 


24.23 


0.38 


1.03 


20.38 


1.75 


4 


X 900 


Q2343-BX436 


23:46:09.06 


12:47:56.00 


2.3277 


23.07 


0.12 


0.47 


21.04 


0.73 


4 


X 900 


Q2343-BX442 


23:46:19.36 


12:47:59.69 


2.1760 


24.48 


0.40 


1.14 


19.85 


2.36 


5 


X 900 


Q2343-BX461 


23:46:32.96 


12:48:08.15 


2.5662 


24.40 


0.44 


0.90 


21.67 


1.99 


4 


X 900 


Q2343-BX474 


23:46:32.88 


12:48:14.08 


2.2257 


24.42 


0.31 


1.15 


20.56 


1.73 


4 


X 900 


Q2343-BX480 


23:46:21.90 


12:48:15.61 


2.2313 


23.77 


0.29 


1.03 


20.44 


1.92 


4 


X 900 


Q2343-BX493 


23:46:14.46 


12:48:21.64 


2.3396 


23.63 


0.28 


0.78 


21.65 


0.75 


2 


X 900 


Q2343-BX513 


23:46:11.13 


12:48:32.14 


2.1092^ 


23.93 


0.20 


0.41 


20.10 


1.87 


4 


X 900 


Q2343-BX529 


23:46:09.72 


12:48:40.33 


2.1129 


24.42 


0.20 


0.91 


21.41 


1.52 


2 


X 900 


Q2343-BX537 


23:46:25.55 


12:48:44.54 


2.3396 


24.44 


0.23 


0.72 


21.43 


1.65 


4 


X 900 


Q2343-BX587 


23:46:29.18 


12:49:03.34 


2.2430 


23.47 


0.32 


1.12 


20.12 


1.82 


3 


X 900 


Q2343-BX599 


23:46:13.85 


12:49:11.31 


2.0116 


23.50 


0.10 


0.81 


20.40 


1.19 


4 


X 900 


Q2343-BX601 


23:46:20.40 


12:49:12.91 


2.3769 


23.48 


0.22 


0.75 


20.55 


1.50 


4 


X 900 


Q2343-BX610 


23:46:09.43 


12:49:19.21 


2.2094 


23.58 


0.34 


0.75 


19.21 


2.24 


4 


X 900 


Q2343-BX660 


23:46:29.43 


12:49:45.54 


2.1735 


24.36 


-0.09 


0.45 


20.98 


2.26 


2 


X 900 


Q2343-MD59 


23:46:26.90 


12:47:39.87 


2.0116 


24.99 


0.20 


1.47 


20.14 


2.59 


4 


X 900 


Q2343-MD62 


23:46:27.23 


12:47:43.48 


2.1752 


25.29 


0.21 


1.23 


21.45 


2.23 


4 


X 900 


Q2343-MD80 


23:46:10.79 


12:48:33.24 


2.0138 


24.81 


0.09 


1.33 


21.38 


1.15 


4 


X 900 


Q2346-BX120 


23:48:26.30 


00:20:33.16 


2.2664 


25.08 


0.02 


0.91 








4 


X 900 


Q2346-BX220 


23:48:46.10 


00:22:20.95 


1.9677 


23.57 


0.29 


0.90 


20.82 






4 


X 900 


Q2346-BX244 


23:48:09.61 


00:22:36.18 


1.6465 


24.54 


0.41 


1.16 








4 


X 900 


Q2346-BX404 


23:48:21.40 


00:24:43.07 


2.0282 


23.39 


0.18 


0.41 


20.05 






5 


X 900 


Q2346-BX405 


23:48:21.22 


00:24:45.46 


2.0300 


23.36 


0.08 


0.60 


20.27 






5 


X 900 


Q2346-BX416 


23:48:18.21 


00:24:55.30 


2.2404 


23.49 


0.40 


0.75 


20.30 






3 


X 900 


Q2346-BX482 


23:48:12.97 


00:25:46.34 


2.2569 


23.32 


0.22 


0.90 








4 


X 900 


SSA22a-MD41<: 


22:17:39.97 


00:17:11.04 


2.1713 


23.31 


0.19 


1.31 








15 


X 720 


West-BM115 


14:17:37.57 


52:27:05.42 


1.6065 


23.41 


0.28 


0.36 








10 


X 900 


West-BX600 


14:17:15.55 


52:36:15.64 


2.1607 


23.94 


0.10 


0.46 








5 


X 900 



^ Un, G and Ti magnitudes are AB. 

^ J and Ks magnitudes are Vega. 

'^ Observed with the ISAAC spectrograph on the VLT; previously discussed bv lErb et alj 120031) . 

'^ Q2343-BX513 was observed a second time with a different position angle, and yielded an Ha redshift Zfja 



: 2.1079. 



TABLE 2 
Near-IR Imaging 



Field 


Band 


Exposure time 

(hrs) 


Depth'' 


GOODS-N 


J 


13.0 


24.1 




Ks 


10.3 


22.6 


Q1623 


J 


9.8 


23.8 




Ks 


11.2 


22.3 


Q1700 


J 


10.7 


24.0 




Ks 


11.0 


22.2 


Q2343 


J 


10.7 


24.0 




Ks 


12.1 


22.3 


Q2346 


Ks 


2.6 


21.2 



Approximate 3(7 image depth, in Vega magnitudes. 
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TABLE 3 
Results of SED Fitting 



Object 


r 


E{B - V] 


)a Ag^b 


SFR'^ 


M.d 




(Myr) 




(Myr) 


(Mq yr-i) 


(lOio Mq) 


HDF-BMZ1156 


500 


0.000 


1900 


4 


9.6 


HDF-BX1277 


const 


0.095 


321 


15 


0.5 


HDF-BX1303 


const 


0.100 


1139 


7 


0.8 


HDF-BX1311 


const 


0.105 


255 


34 


0.9 


HDF-BX1322 


const 


0.085 


360 


19 


0.7 


HDF-BX1332 


const 


0.290 


15 


159 


0.2 


HDF-BX1368 


const 


0.160 


404 


37 


1.5 


HDF-BX1376 


const 


0.070 


227 


9 


0.2 


HDF-BX1388 


const 


0.265 


3000 


23 


7.0 


HDF-BX1397 


const 


0.150 


1015 


17 


1.7 


HDF-BX1409 


const 


0.290 


1015 


27 


2.7 


HDF-BX1439 


const 


0.175 


2100 


21 


4.5 


HDF-BX1479 


const 


0.110 


806 


13 


1.0 


HDF-BX1564 


100 


0.065 


360 


9 


3.3 


HDF-BX1567 


50 


0.050 


227 


5 


2.3 


HDF-BX305 


const 


0.285 


571 


53 


3.0 


Q1623-BX376 


const 


0.175 


64 


80 


0.5 


Q1623-BX428 


50 


0.000 


255 


1 


1.1 


Q1623-BX429 


const 


0.120 


227 


23 


0.5 


Q1623-BX432 


const 


0.060 


1278 


6 


0.8 


Q1623-BX447 


500 


0.050 


1434 


5 


4.4 


Q1623-BX449 


const 


0.110 


2600 


9 


2.4 


Q1623-BX452 


const 


0.195 


3000 


14 


4.3 


Q1623-BX453 


const 


0.275 


454 


107 


4.9 


Q1623-BX455 


const 


0.265 


15 


58 


0.09 


Q1623-BX458 


const 


0.165 


571 


55 


3.1 


Q1623-BX472 


const 


0.130 


3000 


11 


3.2 


Q1623-BX502 


const 


0.220 


6 


72 


0.04 


Q1623-BX511 


const 


0.235 


571 


13 


0.8 


Q1623-BX513 


const 


0.145 


454 


46 


2.1 


Q1623-BX516 


const 


0.145 


1800 


28 


5.1 


Q1623-BX522 


const 


0.180 


2600 


24 


6.2 


Q1623-BX528 


const 


0.175 


2750 


44 


12.2 


Q1623-BX543 


const 


0.305 


8 


528 


0.4 


Q1623-BX586 


const 


0.195 


1434 


17 


2.5 


Q1623-BX599 


const 


0.125 


1900 


35 


6.7 


Q1623-BX663 


1000 


0.135 


2000 


21 


13.2 


Q1623-MD107 


const 


0.060 


1015 


4 


0.4 


Q1623-MD66 


const 


0.235 


905 


43 


3.9 


Q1700-BX490 


const 


0.285 


10 


448 


0.4 


Q1700-BX505 


const 


0.270 


1800 


20 


3.6 


Q1700-BX523 


const 


0.260 


255 


42 


1.1 


Q1700-BX530 


50 


0.045 


203 


6 


1.8 


Q1700-BX536 


50 


0.115 


180 


15 


2.8 


Q1700-BX561 


500 


0.130 


1609 


10 


11.5 


Q1700-BX581 


const 


0.215 


35 


70 


0.2 


Q1700-BX681 


const 


0.315 


10 


628 


0.6 


Q1700-BX691 


1000 


0.125 


2750 


5 


7.6 


Q1700-BX717 


const 


0.090 


509 


8 


0.4 


Q1700-BX759 


const 


0.230 


640 


37 


2.4 


Q1700-BX794 


const 


0.130 


454 


25 


1.1 


Q1700-BX917 


200 


0.040 


806 


4 


4.0 


Q1700-MD103 


const 


0.305 


1015 


65 


6.6 


Q1700-MD109 


const 


0.175 


2200 


8 


1.7 


Q1700-MD154 


const 


0.335 


128 


347 


4.4 


Q1700-MD174 


1000 


0.195 


2400 


24 


23.6 


Q1700-MD69 


2000 


0.275 


2750 


31 


18.6 


Q1700-MD94 


const 


0.500 


719 


213 


15.3 


Q2343-BM133 


const 


0.115 


143 


35 


0.5 


Q2343-BX163 


const 


0.050 


1434 


9 


1.3 


Q2343-BX169 


const 


0.125 


203 


46 


0.9 


Q2343-BX182 


const 


0.100 


180 


23 


0.4 


Q2343-BX236 


const 


0.085 


1680 


13 


2.1 


Q2343-BX336 


const 


0.210 


321 


58 


1.9 


Q2343-BX341 


const 


0.210 


102 


50 


0.5 


Q2343-BX378 


const 


0.165 


255 


11 


0.3 


Q2343-BX389 


const 


0.250 


2750 


22 


6.1 


Q2343-BX390 


const 


0.150 


404 


17 


0.7 


Q2343-BX391 


const 


0.195 


64 


25 


0.2 


Q2343-BX418 


const 


0.035 


81 


12 


0.1 


Q2343-BX429 


const 


0.185 


321 


12 


0.4 


Q2343-BX435 


const 


0.225 


1434 


30 


4.4 


Q2343-BX436 


const 


0.070 


321 


33 


1.1 


Q2343-BX442 


2000 


0.225 


2750 


25 


14.7 
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Object 


r 


E{B - V)^ 


Ageb 


SFR'^ 


M*d 




(Myr) 




(Myr) 


(M0 yr-i) 


(lOio Mq) 


Q2343-BX461 


const 


0.250 


15 


86 


0.1 


Q2343-BX474 


const 


0.215 


2750 


26 


7.2 


Q2343-BX480 


const 


0.165 


905 


33 


3.0 


Q2343-BX493 


const 


0.255 


6 


220 


0.1 


Q2343-BX513 


const 


0.135 


3000 


20 


5.9 


Q2343-BX529 


const 


0.145 


404 


14 


0.6 


Q2343-BX537 


const 


0.130 


571 


15 


0.8 


Q2343-BX587 


const 


0.180 


719 


49 


3.5 


Q2343-BX599 


const 


0.100 


1609 


21 


3.3 


Q2343-BX601 


const 


0.125 


640 


36 


2.3 


Q2343-BX610 


1000 


0.155 


2100 


32 


23.2 


Q2343-BX660 


const 


0.010 


2750 


5 


1.4 


Q2343-MD59 


2000 


0.200 


3000 


11 


7.6 


Q2343-MD62 


const 


0.150 


2750 


7 


1.9 


Q2343-MD80 


50 


0.020 


255 


1 


0.6 


Q2346-BX220 


50 


0.055 


227 


4 


1.7 


Q2346-BX404 


const 


0.095 


1800 


22 


4.0 


Q2346-BX405 


100 


0.010 


321 


7 


1.6 


Q2346-BX416 


const 


0.195 


454 


55 


2.5 



'^ Typical uncertainty {crE(B-V)/l^{B ^ ^)) 
^ Typical uncertainty (crAge/Age) =0.5. 
'^ Typical uncertainty ((tsfr/SFR) = 0.6. 
"^ Typical uncertainty \aMjMi,) = 0.4. 



0.7. 











TABLE 4 
Kinematics 










Object 


ZHa 


2^abs^ 




File,'-' 


(kpc) 


(km s-i ) (km 


s-i) (lOi" Mq) 


(lOlO Mq) 


CDFb-BN88 
HDF-BX1055 


2.2615 
2.4899 


2.4865 


2.4959 


2.6 ±0.2 
2.6 ±0.6 


7.4 
3.8 


96l^» 
< 59 


5.3 
< 1.0 






HDF-BX1084 
HDF-BX1085 


2.4403 
2.2407 


2.4392 
2.2381 




7.3 ±0.2 
1.1±0.1 


4.8 
2.4 


1021^ 


5.1 






HDF-BX1086 


2.4435 






1.8 ±0.2 


< 1.9 


95+20 


< 1.6 






HDF-BX1277 
HDF-BX1303 


2.2713 
2.3003 


2.2686 
2.3024 


2.3051 


5.3 ±0.2 
2.6 ±0.5 


5.0 

3.4 


63+?o 


1.3 


0.5 
0.8 


HDF-BX1311 


2.4843 


2.4804 


2.4890 


8.0 ±0.4 


5.5 


88+9 • 


3.2 


0.9 


HDF-BX1322 


2.4443 


2.4401 


2.4491 


2.0 ±0.2 


4.6 


<47 


< 1.1 


0.7 


HDF-BX1332 


2.2136 


2.2113 




4.4 ±0.3 


5.5 


54+^^ 47 ±14 1.2 


0.2 


HDF-BX1368 


2.4407 


2.4380 


2.4455 


8.8 ±0.4 


7.4 


139+« 


6.6 


1.5 


HDF-BX1376 


2.4294 


2.4266 


2.4338 


2.2 ±0.2 


2.9 


96li 


2.2 


0.2 


HDF-BX1388 


2.0317 


2.0305 




5.8 ±0.5 


7.9 


140±- 


11.7 


7.0 


HDF-BX1397 


2.1328 


2.1322 




5.3 ±0.5 


8.6 


125+g 109 


±19 10.3 


1.7 


HDF-BX1409 


2.2452 


2.2433 




8.5 ±0.6 


7.0 


158tl« 


13.3 


2.7 


HDF-BX1439 


2.1865 


2.1854 


2.1913 


8.8 ±0.3 


8.2 


120+1 


8.6 


4.5 


HDF-BX1479 


2.3745 


2.3726 


2.3823 


2.5 ±0.2 


4.3 


mtf. 


0.9 


1.0 


HDF-BX1564 
HDF-BX1567 


2.2225 
2.2256 


2.2219 
2.2257 




8.6 ±0.7 
4.0 ±0.6 


15.4 
2.6 


99tl^ 
< 62 


11.7 
<1.8 


3.3 
2.3 


HDF-BX305 


2.4839 


2.4825 




4.2 ±0.4 


3.8 


140lg 


5.6 


3.0 


HDF-BMZ1156' 


2.2151 






5.4 ±0.4 


<1.9 


1961^4 


<7.1 


9.6 


Q0201-B13 


2.1663 






2.4 ±0.1 


7.4 


62+ig 


2.2 






Q1307-BM1163 


1.4105 


1.4080 




28.7 ±1.2 


6.7 


12btl 


8.2 






Q1623-BX151' 


2.4393 






3.5 ±0.6 


8.4 


98t!^ 


6.3 






Q1623-BX214 


2.4700 


2.4674 




5.3 ±0.4 


4.1 


55tl^ 


1.0 






Q1623-BX215 
Q1623-BX252 


2.1814 
2.3367 


2.1819 




4.8 ±0.3 
1.1 ±0.3 


5.0 
<1.9 


70+15 


2.0 






Q1623-BX274 


2.4100 


2.4081 


2.4130 


9.5 ±0.3 


4.8 


121+1 ■ 


5.5 






Q1623-BX344 


2.4224 






17.1 ±0.8 


4.1 


92i?o 


2.7 






Q1623-BX366 


2.4204 


2.4169 




7.8 ±1.5 


8.4 


W3tll 


7.0 






Q1623-BX376 
Q1623-BX428 


2.4085 
2.0538 


2.4061 
2.0514 


2.4153 
2.0594 


5.3 ±0.7 
2.7 ±0.5 


7.0 

1.4 


26ltlt 


36.6 


0.5 
1.1 


Q1623-BX429 


2.0160 


2.0142 




5.1 ±0.5 


6.7 


571- : 


1.7 


0.5 


Q1623-BX432 


2.1817 






5.4 ±0.3 


4.6 


5411^6 


1.0 


0.8 


Q1623-BX447 


2.1481 


2.1478 




5.6 ±0.3 


5.3 


l7Atll 160 ±22 12.5 


4.4 
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Object 


^Ha 


^abs'" 




i^Ha" 


(kpc) 


a" 
(km s-i ) 


(kms-1) (IQlO M0) 


(lOio Mq) 


Q1623-BX449 


2.4188 






3.5 ±0.9 


< 1.9 


<72 


<2.5 


2.4 


Q1623-BX452 


2.0595 


2.0595 




4.4 ±0.4 


9.6 


1291^° 


12.9 


4.3 


Q1623-BX453 


2.1816 


2.1724 


2.1838 


13.8 ±0.2 


4.1 


61±1 


1.2 


4.9 


Q1623-BX455 


2.4074 


2.4066 




18.8 ±1.1 


5.0 


isrtll 


13.6 


0.09 


Q1623-BX458 


2.4194 


2.4174 




4.3 ±0.5 


8.4 


160+^!^: 


116 ±18 16.8 


3.1 


Q1623-BX472 


2.1142 


2.1144 




3.9 ±0.2 


6.0 


noil" 


110 ±12 5.7 


3.2 


Q1623-BX502 


2.1558 


2.1549 


2.1600 


13.2 ±0.4 


6.2 


75ti 


2.7 


0.04 


Q1623-BX511 


2.2421 






3.4 ±0.3 


6.2 


152t^^ 


80 ±18 11.3 


0.8 


Q1623-BX513 


2.2473 


2.2469 


2.2525 


3.3 ±0.3 


3.4 






2.1 


Q1623-BX516 


2.4236 


2.4217 




5.2 ±1.0 


2.6 


11A+15 


2.7 


5.1 


Q1623-BX522 


2.4757 


2.4742 




2.8 ±0.3 


6.2 


<46 


< 1.0 


6.2 


Q1623-BX528 


2.2682 


2.2683 




7.7 ±0.4 


7.9 


142±il 


76 ±12 12.0 


12.2 


Q1623-BX543 


2.5211 


2.5196 




8.6 ±0.7 


7.4 


148t^^ 


12.6 


0.4 


Q1623-BX586 


2.1045 






5.1 ±0.4 


5.3 


124ll« 


6.2 


2.5 


Q1623-BX599 


2.3304 


2.3289 


2.3402 


18.1 ±0.4 


5.8 


162+1 


11.6 


6.7 


Q1623-BX663' 


2.4333 


2.4296 


2.4353 


8.2 ±0.3 


7.0 


I32tll 


9.3 


13.2 


Q1623-MD107 


2.5373 






3.7 ±0.4 


3.6 


<43 


<0.9 


0.4 


Q1623-MD66 


2.1075 


2.1057 




19.7 ±0.3 


5.3 


120+^ 


5.9 


3.9 


Q1700-BX490 


2.3960 


2.3969 


2.4043 


17.7 ±0.6 


5.5 


110t« 


5.1 


0.4 


Q1700-BX505 


2.3089 






3.6 ±0.3 


7.4 


120+2; 


8.3 


3.6 


Q1700-BX523 


2.4756 






4.7 ±0.5 


7.7 


130t^^ 


9.9 


1.1 


Q1700-BX530 


1.9429 


1.9411 




12.2 ±0.7 


6.2 


< 37 


<0.6 


1.8 


Q1700-BX536 


1.9780 






11.3 ±0.7 


7.9 


89tlt 


4.7 


2.8 


Q1700-BX561 


2.4332 


2.4277 




1.9 ±0.6 


< 1.9 






11.5 


Q1700-BX581 


2.4022 


2.3984 




4.0 ±0.7 


4.3 






0.2 


Q1700-BX681 


1.7396 


1.7398 


1.7467 


6.3 ±0.2 


8.6 






0.6 


Q1700-BX691 


2.1895 






7.7 ±0.3 


6.7 


170+11 


220 ±14 15.0 


7.6 


Q1700-BX717 


2.4353 




2.4376 


3.8 ±0.4 


4.8 


<47 


<1.1 


0.4 


Q1700-BX759 


2.4213 






1.3 ±0.5 


2.4 






2.4 


Q1700-BX794 


2.2473 






6.8 ±0.4 


4.8 


80+11 


2.4 


1.1 


Q1700-BX917 


2.3069 


2.3027 




7.4 ±0.5 


10.6 


9911^ 


8.2 


4.0 


Q1700-MD69 


2.2883 


2.288 




8.2 ±0.9 


9.4 


155111 


17.2 


18.6 


Q1700-MD94' 


2.3362 






2.8 ±0.3 


9.6 


730tll 


382 


15.3 


Q1700-MD103 


2.3148 






4.1 ±0.8 


8.4 


75+21 
'0_24 


100 ±19 3.7 


6.6 


Q1700-MD109 


2.2942 






8.9 ±0.7 


4.1 


93+^^ 


2.7 


1.7 


Q1700-MD154 


2.6291 






7.5 ±0.5 


8.2 


57tt^ 


2.1 


4.4 


Q1700-MD174' 


2.3423 






12.9 ±1.4 


3.6 


4441^0 


56.5 


23.6 


Q2343-BM133 


1.4774 


1.4769 




28.7 ±0.8 


7.7 


t>0_6 


1.8 


0.5 


Q2343-BM181 


1.4951 


1.4952 




3.4 ±0.5 


7.7 


391JS 


45 ±18 1.0 




Q2343-BX163 


2.1213 






2.2 ±0.4 


4.1 






1.3 


Q2343-BX169 


2.2094 


2.2105 


2.2173 


4.7 ±0.3 


3.1 






0.9 


Q2343-BX182 


2.2879 


2.2857 


2.2909 


2.4 ±0.3 


2.6 






0.4 


Q2343-BX236 


2.4348 


2.4304 


2.4372 


3.1 ±0.6 


4.8 


148+1" 


8.2 


2.1 


Q2343-BX336 


2.5439 


2.5448 


2.5516 


4.3 ±0.6 


5.0 


I40I34 


7.6 


1.9 


Q2343-BX341 


2.5749 


2.5715 




4.0 ±0.6 


2.4 






0.5 


Q2343-BX378 


2.0441 






4.5 ±1.0 


7.4 






0.3 


Q2343-BX389 


2.1716 


2.1722 




12.0 ±0.4 


10.1 


lllll 


71 ±10 9.6 


6.1 


Q2343-BX390 


2.2313 


2.2290 




4.9 ±0.5 


10.8 


78tll 


5.1 


0.7 


Q2343-BX391 


2.1740 


2.1714 




4.2 ±0.2 


6.0 






0.2 


Q2343-BX418 


2.3052 


2.3030 


2.3084 


8.0 ±0.2 


3.6 


66+1 


1.2 


0.1 


Q2343-BX429 


2.1751 






4.8 ±0.3 


9.4 


01-20 


1.9 


0.4 


Q2343-BX435 


2.1119 


2.1088 


2.1153 


8.1 ±0.4 


8.9 


60111 


2.4 


4.4 


Q2343-BX436 


2.3277 


2.3253 


2.3315 


7.2 ±0.4 


9.1 


63+1° 


2.8 


1.1 


Q2343-BX442 


2.1760 






7.1 ±0.3 


10.6 


132+^ 


13.5 


14.7 


Q2343-BX461 


2.5662 


2.5649 


2.5759 


7.0 ±0.7 


7.7 


139li 


11.5 


0.1 


Q2343-BX474 


2.2257 


2.2263 




5.0 ±0.3 


7.0 


S/1 + 12 
84_i2 


3.8 


7.2 


Q2343-BX480 


2.2313 


2.2297 


2.2352 


3.0 ±0.2 


3.4 


< 37 


<0.7 


3.0 


Q2343-BX493 


2.3396 


2.3375 


2.3447 


5.3 ±0.9 


6.5 


155+il 


12.2 


0.1 


Q2343-BX513J 


2.1092 


2.1090 


2.114 


10.1 ±0.4 


4.1 


150tl 


7.0 


5.9 


Q2343-BX529 


2.1129 


2.1116 


2.1190 


3.5 ±0.5 


<1.9 






0.6 


Q2343-BX537 


2.3396 






5.2 ±0.3 


5.3 






0.8 


Q2343-BX587 


2.2430 


2.2382 




5.5 ±0.4 


7.7 


94+- 


5.2 


3.5 


Q2343-BX599 


2.0116 


2.0112 




4.5 ±0.4 


7.2 


771^1 


3.3 


3.3 
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Object 



ZHo 



^abs 



^hya 



i^Ha" 



(kpc) (km s-i ) (km s-i ) (10^° M©) (lO"' Mq) 



Q2343-BX601 


2.3769 


2.3745 


2.3823 


7.4 ±0.4 


5.3 


10511^ 




4.5 


2.3 


Q2343-BX610 


2.2094 


2.2083 


2.2129 


8.1 ±0.4 


8.4 


96t?o 




6.1 


23.2 


Q2343-BX660 


2.1735 


2.1709 


2.1771 


9.4 ±0.4 


6.2 


< 40 




<0.8 


1.4 


Q2343-MD59 


2.0116 


2.0107 




2.9 ±0.4 


5.5 








7.6 


Q2343-MD62 


2.1752 


2.1740 




2.3 ±0.3 


6.7 


79+27 
'^-30 




3.3 


1.9 


Q2343-MD80 


2.0138 


2.0116 




3.2 ±0.3 


2.6 


74+1' 




1.1 


0.6 


Q2346-BX120 


2.2664 






5.3 ±0.3 


6.2 


621- 


40 ±10 


1.9 




Q2346-BX220 


1.9677 


1.9664 




10.3 ±0.6 


6.5 


143111 




10.4 


1.7 


Q2346-BX244 


1.6465 


1.6462 


1.6516 


5.4 ±0.8 


7.7 


42l^^ 




1.0 




Q2346-BX404 


2.0282 


2.0270 


2.0348 


13.9 ±0.3 


3.6 


102l| 




3.0 


4.0 


Q2346-BX405 


2.0300 


2.0298 


2.0358 


14.0 ±0.2 


5.8 


5011 




1.1 


1.6 


Q2346-BX416 


2.2404 


2.2407 




12.1 ±0.7 


4.3 


12611^ 




5.3 


2.5 


Q2346-BX482 


2.2569 


2.2575 




11.2 ±0.3 


9.8 


133+1 




13.7 




SSA22a-MD41 


2.1713 






7.9 ±0.1 


9.6 


107l« 


150 ± 16 


8.2 




West-BM115 


1.6065 


1.6060 




5.9 ±0.4 


5.5 


12811^ 




6.9 




West-BX600 


2.1607 






6.3 ±0.4 


8.2 


181+1« 


210 ±13 


21.5 





'^ Vacuum rcdshift of the UV insterstellar absorption lines. We give a value only when the S/N of the spectrum is sufficient for a precise 
measurement. 

^ Vacuum rcdshift of the Lya emission line, when present. 

'^ Observed flux of Ha emission line, in units of 10~^^ erg s~^ cm~2 

'^ Approximate spatial extent of the Ha emission (FWHM), after subtraction of the seeing in quadrature. 

° Velocity dispersion of the Ho emission line. 
For tilted emission lines, the velocity shear (umax — fmin)/2, where Dmax and Dmin ^rc with respect to the systemic rcdshift. 

s Dynamical mass M^y^ = 3A<j^riia/G. 

^ Stellar mass, from SED modeling. 

' AGN 

J Q2343-BX513 was observed twice with NIRSPEC, with position angles differing by 9°. The first observation yielded ZHa = 2.1079 
and (T = 58 km s^^ , and the second ^jja = 2.1092 and a = 150 km s^^ . It also has two Lya emission rcdshifts, ZLya = 2.106 and 

2Ly« = 2.114. 



